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PRESENTATION 


Itis with great enthusiasm and dedication that | present to you, reader, the 
new collection “Engineering in perspective science, technology and innovation”. 
A collection of scientific works that addresses current and highlighted topics in 
the field of Engineering. 

This consists of four chapters, covering various subjects carried out in 
various teaching and research institutions abroad in Brazil. Relevant issues such 
as the potential for reusing slate waste in the ceramics industry; the Integration of 
E-learning in the metaverse for the development of teaching and training in higher 
education; Thermal study of a coated oxygen for the production of hydrogen, as 
well as thermal-hydraulic advances in the production of this material. 

A new, current and contemporary collection, aiming to provide learning 
opportunities for readers in both the academic and professional fields. 
Furthermore, this work seeks to encourage the dissemination of new work and 
highlights the importance of researchers in publishing their work through reliable 
platforms, such as Atena Editora. 

To the authors, | thank you for your trust and spirit of partnership. 

Good Reading 


Mariana Natale Fiorelli Fabiche 
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ABSTRACT: Throughout the years, mining 
activity has had great importance for 
economic progress, despite its repercussions 
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of generating significant amounts of waste, 
impacting sustainability and environmental 
conservation. The present study aimed to 
investigate the feasibility of using slate waste 
in ceramic processing, specifically through 
slip-casting techniques. By conducting 
detailed experimental procedures and 
characterization tests, the study seeks to 
provide insights into the suitability of slate 
waste as a substitute for conventional raw 
materials. Additionally, the study aimed 
to evaluate the technical properties and 
performance of ceramic pieces produced 
using slate waste, thereby addressing 
concerns about the quality and applicability 
of such materials in practical applications. 
Slate rejects collected from quarries were 
characterized and ceramic pieces were 
produced by slip casting. The cast green 
and sintered samples were characterized 
by density, pore size distribution by 
mercury porosity, surface area evaluation, 
X-ray diffraction techniques, and scanning 
electron microscopy. The X-ray results, 
when confronted with thermal analyses, 
have shown major structural transitions 
at approximately 550ºC and 850ºC. The 
studies of porosity and density of pieces, 
before and after heat treatment, provided a 
better understanding of the slate sintering 
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process. The results indicate that the density decreases after heat treatment and the pieces 
show a volume reduction due to slate reactions. Compacts formed by slip-casting different 
suspensions showed an average pore size of approximately 0.72 mm, suggesting that the 
structural changes that occurred in the slate greatly influenced the porosities and densities of 
the resulting materials. The results allow us to infer that slate powder from rejects has a great 
potential for being recovered, recycled, and applied in ceramic processing, producing various 
ceramic products. These products encompass a wide range of applications, including but not 
limited to sanitary ware, decorative pieces, and architectural elements. 

KEYWORDS: Slate waste, Material reuse, Recycling, Ceramic processing, Sustainability. 


INTRODUCTION 


Different society segments have been devoting attention to slate waste treatment 
techniques and recycling, especially universities, tech institutes, and industries, focusing 
on the massive amount of waste produced in its extraction sites. These rejects are almost 
always incorrectly disposed of by the respective companies, contributing to environmental 
degradation. 

Brazil is accounted for its varied geodiversity, supporting its worldwide importance 
as an exporter of ornamental rocks. Truly, the country is the third largest slate producer 
and exporter, with the state of Minas Gerais accounting for about 90% of this production 
(approximately 0.4 million tons/year) and almost all Brazilian exports. Slate derivates 
from Minas Gerais unfold into approximately 18 million square meters of slabs, tiles, pool 
countertops, roofs, and other products, according to Abirochas (2018). The natural slate 
extraction, cutting, and polishing process is extremely wasteful and its reject stream from 
quarrying is not properly managed in the present day. According to Oti, Kinuthia and Bai 
(2010) finished pieces derived from slate have shape and size requirements that make 
it unfeasible to use all extracted raw material, so only around 10% end up entering the 
production process, and even less reach the final phases. Frias et al. (2014) report that 
slate used in civil construction generates between 75% and 90% in weight of waste, 
and one possible way of reusing it is through thermal activation and manufacturing new 
cementitious materials. The large amounts of engendered slate waste happen because 
of many variables: poor extractive and processing activity technological profile; research 
deficiency on qualifying new mining fronts; lack of boulder-type slate exploitation; resistance 
to directing sales to the foreign market as long as the domestic one, which would be more 
profitable allied to public-private partnerships to promote less bureaucratical processes, to 
name a few. All of those result in an excessive quantity of rejects — small pieces of irregularly 
shaped or liquid-effluent rocks, and also mud from cutting and polishing processes (Palhares, 
Santos and Hunter, 2016). 

The waste produced annually by Brazilian slate industries is nearly 0.3 million tons 
(Oti, Kinuthia and Bai, 2010). Itis composed mainly of water, lubricants, crushed rock, and 
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slate fragments. The production of alternative materials with the slate waste brought as a 
constituent to the rock manufacturing industries may reduce or even eliminate pollution in 
the extraction areas, aside from promoting new job opportunities and income — essential 
to the country's progress and development. Many studies have reported the production of 
new materials with slate waste through traditional ceramic processes, such as pressing, 
extrusion, slip casting, and others. In Palhares et al, (1997) work, there were several slate 
bricks produced by pressing and heat-treating the samples up to 1000ºC. Their compressive 
strength tests showed that bricks have great potential for use in civil construction. Catarino 
et al. (2003) used residue powders from cutting and machining operations in slate quarries 
to produce flooring and roofing tiles by a press (at 40 MPa) and sintering route (up to 
1150 ºC). The authors showed that the properties of the products obtained were suitable 
for flooring tiles. Cambronero, Ruiz-Roman and Ruiz-Pietro (2005) obtained foams (porous 
ceramic) without any foaming agent addition and, according to them, the resulting material 
reaches adequate strength after the transformations suffered by the slate components 
during heating. Silva and Peres (2006) investigated the pozzolanicity of slate wastes 
and confirmed such activity with calcinated slate samples but not with natural ones. They 
proposed the use of slate wastes in the cement industry as mineral addition in partial 
substitution to Portland cement clinker. Snelson et al. (2010) reviewed the potential of using 
slate waste and suggested the great potential for using slate waste in many situations, 
including engineering applications despite reported transportation difficulties. Generally, the 
use of slate wastes is feasible and could minimize the environmental impact of their disposal 
but in countries with a large territorial area such as Brazil, without an adequate railway 
network, this transport can increase costs and make their use impracticable for ceramic 
industries. Zhen, et al, (2018) processed the slate wastes to form ceramsite slate, a kind of 
artificial foam, similar to that obtained by Cambronero, Ruiz-Roma and Ruiz-Pietro (2005). 
The ultra-light ceramsite could be directly used as building aggregate since its properties 
were shown to be superior to other ceramsites. At present, mercury intrusion porosimetry 
(MIP) is still today the reference experimental method to obtain the pore size distribution 
of dry samples of porous media (Cieszko, Kempinski and Czerwinski, 2019), despite its 
environmental health and safety concerns. Some substantial benefits of this extensively 
used technique include the relatively fast measurements, the possibility of covering a large 
range of pore sizes, i.e., ffom 10º up to 10º m (Zhang, et al, (2016), and the richness of 
accessible databases. 

Inthe present work, the utilization of slate waste in ceramic processing was examined, 
emphasizing its feasibility as an alternative resource and its potential contributions to 
sustainability and resource conservation. This initiative aimed to enhance economic efficiency 
while concurrently addressing the issue of overexploitation of clay quarries. The pore 
structure of ceramic pieces manufactured via slip casting was evaluated. Characterization 
tests, utilizing Mercury Intrusion Porosimetry (MIP) techniques, were conducted before and 
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after burning to investigate morphology modifications and pore structure changes during 


thermal treatment. 


EXPERIMENTAL PROCEDURE 


The slate powder used in this work was kindly donated by the Company Micapel 
Slate, located in the city of Pompeu, state of Minas Gerais, and was produced via block 
extraction for the production of slabs, floors, and roofs. 

It was pre-treated before characterization to remove impurities such as organic 
waste and lubrification contamination. The steps of this treatment were: preparing for 
powder suspension in distilled water, wet sieving with the 4400 sieve; 24-hour decanting; 
siphonation (water removal), and oven drying at 120ºC for 24 hours. 

Powder composition and crystallinity were determined with a diffractometer XRD 
Philips - X-Ray Generator, No. DY 1308, type PW 1730/10 with Cu-K, rad radiation, 35kV, 
goniometer speed of 0,02º/min in 20 per step, with a 5 seconds per step time count and 
collection from 5º to 80º in 20. For the refinement, the General Structure Analysis System 
Program (Larson and Von Dreele, 2001) was employed with the EXPGUI interface, using the 
Thompson-Cox-Hastings pseudo-Voigt profile function, with background radiation adjusted 
by the Chebyshev polynomial. The mineralogical analysis of the slate powder was obtained 
through refinement by the Rietveld method. The adjustment parameters c2, R, and R,, 
for the analyzed samples were 5.5%, 5.6%, and 6.6%, respectively. The semi-quantitative 
analyses were carried out on a Shimadzu 7000, model EDX-720 equipment. 

SEM images were obtained with a Scanning Electron Microscope LEO/Zeiss 1530 
(LEO Elektronmikroskope GmbH, Germany) joined with energy dispersive spectroscopy 
(EDS). The slate particles were suspended in water and a few drops of these suspensions 
were added to the stub and dried at 60ºC. After that, the samples were uniformly coated 
with gold. 

Thermogravimetric analyses were performed on a Shimadzu TA 50-WSI heat 
analyzer. The analysis was carried out in the temperature range of 20 “C-1150 ºC with a 
heating rate of 10 K min”! under air atmosphere. The slate powder suspensions in distilled 
water were prepared with solids percentages varying from 48%w/v to 78%w/v (suspensions 
1 to 3) and homogenized with a magnetic stirrer for 24 hours. Ammonium polyacrylate 
(pureness degree — 98%) was used as a dispersant (varying from 1%v/v to 2.5%v/v — 
suspensions 4-6) and poly (vinyl alcohol) as a binder (pureness degree — 99.8%, varying 
from 1%v/v to 2.5%v/v — suspensions 7-9); both chemicals were purchased from Aldrich 
and used as received. 

For the manufacture of plaster molds, preforms were used and cut into sizes of about 
2 inches in length. The molds were obtained with plaster to water ratio of 70 (100 g gypsum 
to 70 mL of water), and the mixture was stirred for 1 minute and poured into the pre-mold 
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to form the wish piece structure. After curing the mold, it was dried at room temperature for 
15 days before use. 

After casting, the pieces were dried at room temperature for 24 hours, followed by 
oven heating at 100'C and finally heat-treated in an electric furnace at the temperatures 
500ºC, 700ºC, 1000ºC, 1050ºC and 1100ºC with a 10ºC/min heating rate in air. The green 
sample is generally referred to as “cast” since itis measured before heat treatment. The as- 
cast and after-sintering densities were determined by a helium pycnometer (Multipycnometer 
— Quantachrome). The Auto Pore Ill 9420 (Micromeritics) was used to measure the mercury 
penetration porosimetry, for determining pore size distributions. Both low-pressure (0 to 345 
kPa) and high-pressure (345 MPa to 414 MPa) penetration volumes were measured and 
collected automatically. A computer program was used to combine the low-pressure and 
high-pressure data and to convert the volume-pressure data to volume-diameter data. The 
BET method (Autosorb-1, Quantachrome), based on the N, adsorption model (at 77 K), was 
used to evaluate the average surface area. 


RESULTS AND DISCUSSION 


Powder Characterization 


SEM images of the slate powder (Figure 1) show agglomerated particles having 
sizes ranging from 10 to 30 mm. There were different shapes identified among the particles, 
isotropic and laminar, with the appearance of leaves or small plates due to the large variety 
of minerals present in the slate. 

The density values found via helium pycnometer and surface area analysis via 
nitrogen adsorption were, respectively, 2.76 g.cmº* and 12.3 m?, g'. 

EDX analysis of the powder indicated abundant contents of SiO, (61.6%), ALO, 
(16.6%), K,O (10.8%) and Fe,O, (5.9%). These four oxides are the major constituents of 
slate. Catarino et al. (2003) found comparable values using the chemical analysis by FRX 
and reported the presence of 54.68% SiO,, 9.95% Fe,O,, and 23.52% ALO,. According to 
Cambronero, Ruiz-Roman and Ruiz-Pietro (2005) the chemical composition of the slate 
powder consists mainly of SiO, (50-60%) and ALO, (20-25%) and the XRD analysis of 
the powder showed peaks corresponding to the following constituents: quartz, hematite, 
chlorite, and muscovite. 

The XRD analysis of slate powder identified the components: (Q) Quariz - SiIO, 
(29,6%), (C) Chlinoclore (19,2%), (M) muscovite (29,0%), (A) Albite (14,5%), (Ca) Calcite 
(1.8%), and (O) Orthoclase (5,9%) with errors 0,5; 0,4:0,4;0,7:0,7;0,4 respectively. 

Comparing the chemical and mineralogical analyses, it is generally observed that 
silica (SiO,) is the dominant mineral present in quartz, phyllosilicates (micas and chlorites), 
and feldspars. Alumina (Al,O.) is mainly contained in phyllosilicates and, to a lesser extent, 
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in feldspars, the same way as sodium and potassium oxides. Part of the potassium is 
present in muscovite, and the magnesium oxides originate from the clinochlore and the 
magnesium calcite. 

Brazilian slate samples reported elsewhere displayed values similar to those found in 
this work. The values were 26-30% quartz, 32-34% mica, 18-20% chlorite, 12-15% feldspar, 
and 0.5-1% carbonates. Campos, Velasco e Martinez (2004), Chiodi Filho, Rodrigues 
e Arthur (2003) and Frias et al. (2014) also mention the presence of chamosite, quartz, 
muscovite, and feldspar. 

The values obtained for the parameters c2, R,. and Aa refinements all fall within the 
range observed for natural multi-mineral systems, as demonstrated by Hill, Tsambourakis 
and Madesen (1993), Mumme et aí. (1996) and Weidler et al. (1998). 

The thermal analysis showed the main transitions of slate powder. At transition at 
578ºC, it may be assigned to the allotropic transformation of quartz and clinochlore and 
muscovite dihydroxylation. Around 740'C, the muscovite and clinochlore elements lose their 
lattice structure, and around 889ºC, the calcite is decomposed. The mass losses associated 
with the described events were, respectively, 1.5%, 3.2%, 1.6%, and 0.87%. 


Characterization of pieces 


Table 1 shows the bulk density and volume reduction for the pieces obtained. 
Suspensions 1 to 3 confirm that when the solids percentage increase, the green density 
increases, due to the larger number of particulates in suspension, leading to less water 
release during the casting process. Conversely, when the dispersant percentage increased 
(suspensions 4 to 6), the green density decreased. This is explained by the particles 
being further apart, due to the increasing amount of dispersant used. When comparing 
the suspensions 1 to 3 with 4 to 6, the density found is higher in the latter set. In these 
cases, the dispersion is more effective inducing a slower deposition of the particle, which 
renders filling the empty spaces more effective, resulting in higher densities. This suggests 
that the binder could be competing with the dispersant in the surface particle, resulting in 
a less effective dispersion and an attraction between the slate particles, leading to lower 
densities. The use of the binder to obtain “as cast” pieces with adequate handling strength 
was adverse to the bonding process. For suspensions 7 to 9, the density is always higher 
than other suspensions due to the binder, which promotes particle bonding. By increasing 
the binder fraction, the green density decreases. It was also observed density reduction for 
all suspensions subsequently to the heat treatment. 

When assessed, the density values were distant from the natural slate stone 
theoretical estimates, and these deviations could be caused by the thermal expansion 
of clinochlore and muscovite constituents. Along with it, there was an 8-10% weight loss 
caused by outgassing and volatile decomposition. 
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The pore distribution curves for the slate samples obtained through mercury 
porosimetry after casting, across all suspensions, reveal a tendency towards the 
macroporous range with a higher average size, corroborating with the density values. 

The pore size distribution curves before and after heating of the suspensions 
prepared with different solids percentages and with constant dispersant and binder showed 
similar behavior. Before the heat treatment, the curves are not monomodal and the pore 
size range is relatively wide. The highest porosity range for green pieces is around 0.1 to 
imm. As the percentage of solids increased, the average pore diameter decreased from 
approximately 0.72 to 0.45mm, suggesting that the particle cluster was different and more 
efficient when the percentage of solids was increased. These values are consistent with the 
density growth, which in turn occurs along with the increase of solids percentage. 

After the heat treatment, the distribution becomes relatively narrower and shifts to an 
average pore size of approximately 0.72 mm. This may be explained by crystal transitions and 
slate dihydroxylation. This process is not homogeneous, and it occurs in two steps. The first 
stage occurs at about 500ºC and involves transitions related to interlayer water loss, whereas 
the second phase is caused by dihydroxylation and takes place at approximately 700ºC 
(Guggenheim, Chang and Van-Gross, 1987). Furthermore, at high temperatures, changes in 
lattice structure could be taking place so that the tetrahedral sheets become hexagonal and fully 
extended onto their sides, resulting in structure expansion (Hazen, 1997). 

For all samples, the porosity observed varied between 1 and 0.1 mm after the 
heat treatment. This suggests that changes in additive and solids percentages did not 
have an accentuated effect on the pore size distribution of the pieces and the structural 
transformations on the slate are primarily responsible for the bulk residual porosity. The 
thermal treatment reduces the total porosity but not quite the pores' size. These results are 
consistent with some clays that burn up to 1000ºC and show a porosity increase as a result 
of the loss upon ignition, related to the decomposition of clay minerals and organic materials 
(Shu, et al., 2012). 

Usually, higher densities are related to smaller volume reductions; however, in this 
case, the green density is not the only factor influencing the volume shrinkage after the heat 
treatment. The structural changes that occurred in the slate at higher temperatures have a 
significant effect on the final properties of the ceramic pieces obtained through slip casting. 

The structural changes describe the complex thermal transformations of the phase 
minerals of which the slate is formed. Many authors have studied the slate and reported 
its mineral reactions. Sanchez-Soto et al., (2010) reported that chlorite and muscovite 
disappear at around 800ºC and quartz is undragged, making this process exothermic. An 
increase in expansion up to 1020ºC was interpreted by a bloating effect, caused by trapped 
gases in the pores and the release of oxygen gas, being the latter characteristic of raw 
materials containing iron oxides. All reactions promote a piece expansion, reducing the final 
density and increasing porosity. 
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According to Moura and Grade (1977) slate residual pores are due to the thermal 
expansion resulting from gases formation - CO,, SO, and H,O - accumulated in gas 
pockets, gases which are derived from the decomposition of organic matter, hydrous fine- 
grained phyllosilicates, such as illite, chlorite, muscovite, calcite, and hematite. Some of the 
evolved gases could become sealed inside the inorganic mass constituting the fired slate 
piece. The reactions that promote slate expansion take place within the temperature range 
of 400ºC-820ºC, and essentially consist of the aforementioned evolved gases; from 820ºC 
onwards the vitrification phase is initiated and progressively developed, resulting in the 
formation of an essentially alkaline alumino-silicate glass that provides gas retention in the 
most internal parts of black slate pieces. The slate expansion increases up to 1050ºC-1150 
ºC, being the highest when the temperature reaches the maximum temperature value. As 
far as temperature increases and provides gas dilatation, this gas phase presses the glassy 
material that acquires and develops plasticity and elasticity. Gases can migrate up to the 
external surfaces of the slate pieces and escape to the kiln atmosphere (Costa, Almeida e 
Gomes, 2013). 

The density and pore size of green compacts produced through slip casting are 
also influenced by the dispersion state of powder particles in suspension. The development 
of agglomerates may be controlled by adjusting the pH. Meanwhile, the dismantling of 
agglomerated particles may help obtain pieces with higher density and lower porosity. In 
addition, the thermal treatment at higher temperatures may form a liquid phase and fill the 
existing voids in the slate pieces, reducing the initially present porosity in the green fraction. 

Table 1 shows the porosity of the samples before (green samples) and after heat 
treatment. As expected, all samples showed a porosity reduction after heat treatment. The 
remaining porosity is due to the structural modifications suffered by the slate, as discussed 
above. 

The ceramic slate pieces porosities presented a substantial value to all suspensions. 
Suspension 7 offered the best values, which can be associated with a more dispersed 
suspension and lower binder and dispersant concentrations, hence reducing agglomerates. 
In the case of suspensions formed by different minerals, their stability depends on each 
component's physicochemical surface characteristics. Their contribution depends on their 
relative quantity as well as their surface reactivity towards the other elements of the system. 
The minerals present in the slate have very different surface characteristics regarding their 
acidity/basicity concerning water. Therefore, it is inferred that suspension 7 was the most 
stable among those investigated in this work. 

The series of SEM images in Figure 2 exhibits the microscopic evolution of the pieces 
obtained after heat treatment of suspension 7 at temperatures from 500ºC to 1100ºC. Figure 
3 illustrates the corresponding evolution of the diffractograms. 

From the SEM analysis, it was noticed that after burning, the samples acquired a 
smoother surface, making them less porous and stronger than the respective green pieces. 
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At elevated sintering temperatures, the presence of a liquid phase was perceived, but it was 
not enough to fill the pores (Costa, Almeida, and Gomes, 2013). 

From X-ray diffraction results, it was possible to identify important reactions that 
occurred in slate pieces. New phases were identified, such as gehlenite, hercynite, and 
spinels. According to Binda et al. (2020), the presence of these minerals was anticipated, 
due to the presence of calcite in the raw material. The gehlenite is formed by reacting 
calcium oxide with Si and Al from the clay decomposition. At temperatures up to 500ºC, 
the slate powder does not change. Starting at 700ºC, the peaks referring to muscovite 
and clinochlore begin to disappear due to the release of structural water present in these 
minerals and the modifications in their crystalline structures. This observation is confirmed 
using differential thermal analysis with peaks between 700ºC and 800ºC. At temperatures 
above 1000ºC, the presence of the amorphous phase is evidenced by the area under the 
peaks, due to a low-melting liquid phase formation (Palhares et al., 2006). 

The rise in temperature during the thermal process provides a reduction in surface 
area and the number of pores and voids (Table 2). 


CONCLUSIONS 


An experimental study was conducted to investigate the use of slate waste powder 
to produce ceramic pieces via the slip-casting process. It evaluated the influence of residue 
on the density and porosity of the pieces in temperatures until 11000C. 

Employing the slip-casting technique with slate waste as a raw material opens up 
avenues for the production of various ceramic products. These products encompass a wide 
range of applications, including but not limited to sanitary ware, decorative pieces, and 
architectural elements. The utilization of slate waste not only offers economic benefits but 
also contributes to sustainability by reducing the exploitation of clay quarries and facilitating 
the recovery of slate waste. By harnessing this alternative resource, the ceramic industry can 
diversify its product portfolio while minimizing its environmental footprint. The investigation 
underscores the potential of slate waste as a viable raw material in ceramic processing, 
highlighting its significance in promoting resource conservation and advancing sustainable 
practices within the industry. 

The major constituents of slate powder were identified by X-ray diffraction as quartz, 
muscovite, and clinochlore not presenting different components of slate stone. 

The studies of porosity and density of slate ceramic pieces, before and after heat 
treatment, provided a better understanding of the slate sintering process. The results 
indicate that the bulk density decreases after heat treatment and the pieces show a volume 
reduction due to slate reactions. Compacts formed by slip-casting different suspensions 
showed an average pore size of approximately 0.72 mm, suggesting that the structural 
changes that occurred in the slate greatly influenced the porosities and densities of the 
resulting materials. 
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Based on XRD results and thermal analysis of the slate samples and burned 
ceramic pieces, it was possible to identify several important alterations during the burn, 
which may have contributed to the porosity surge. Evaporation of free water, both physically 
and chemically adsorbed (at 250ºC and 320ºC, respectively), dihydroxylation of clays with 
the elimination of structural water (muscovite and clinochlore, 570ºC-580ºC), and calcite 
decomposition with CO, release (800ºC-950ºC) are all events which increase porosity. 

The use of slate as a raw material has great potential due to its physical and 
chemical properties, them being similar to natural clays. lts powder may be used: in 
extrusion processes to produce tiles and bricks; as a filler in various ceramic processing; 
as a concrete aggregate; and, if heat treated, may be employed as a light aggregate, in the 
production of artificial stone materials for many applications. 
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Sample Green bulk : Density after ola Green porosity Porosity after 
number density (g/em?) heat treatment reduction (Y%) (%) + 0,1 heat treatment 
suspension + 0,1 (g/cm?) + 0,1 (1) (%) 0,1 

1 2,1 1,9 17 56,9 46,9 

2 24 2,2 19 56,5 45,7 

3 2,5 2,4 19 55,4 43,5 

4 2,7 2,6 23 59,6 44,6 

5 2,5 2,4 23 63,4 45,4 

6 2,4 2,2 22 64,7 41,1 

7 3,5 2,7 24 50,5 39,8 

8 3,4 2,6 25 50,9 43,4 

9 3,1 2,4 26 52,9 46,5 


Table 1: Green porosity/density of pieces and porosity/density after heat treatment, 


Temp. (ºC) Specific surface BET (m?/g) Total Vol. of Pores (cm*/g) 
100 12,26 0,0378 
500 12,78 0,0365 
700 9,99 0,0353 
1000 4,76 0,0059 
1100 4,22 0,0054 


Table 2: Porosity and surface area of pieces after burning. 
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INTEGRATING E-LEARNING IN METAVERSE 
FOR SUSTAINABLE DEVELOPMENT OF HIGHER 
EDUCATION AND TRAINING 


Nguyen Duc Son 


ABSTRACT: Today, the development of 
information technology in the context of the 
industrial revolution 4.0 has changed the face 
of education rapidly and raised many new 
problems in this context. Along with that, the 
process of globalization has also given birth 
many non-traditional education methods. 
These types of education have transcended 
space, time, and types of boundaries. The 
article mentions the role and importance 
of metaverse in the global trend of higher 
education development through e-learning 
in the process of implementing cross- 
border education in universities all over the 
world. The article uses an interdisciplinary 
approach: education, science, and 
technology... to build a metaversity model in 
metaverse to develop training at universities 
in a sustainable way. 

KEYWORDS:  metaverse, metaversity, 
e-learning, borderless education. 
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INTRODUCTION 

Entering the 21% century, in the 
context of the industrial revolution 4.0, 
applying 
information and communication technology 


higher education has been 
to serve higher education. This is reflected 
in the digital transformation taking place 
in training institutions around the world. 
According to Susan Grajek and Betsy 
Reinitz [13], transfer process at universities 
around the world by 2020 indicates that 13% 
of colleges and universities are engaged in 
digital transformation, 32% are developing 
a digital transformation strategy and 38% 
of other higher education institutions are 
exploring digital transformation, only 17% 
of educational institutions do not have time 
to enter the process. 

Along with that, the process of 
globalization has given birth to many non- 
methods. New 


traditional educational 


generations demand different delivery 
methods and curricular content. Flexibility 
and personalized offerings have become 


a common denominator in all sectors [11]. 
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These types of education have transcended space, time, and types of boundaries. Many non- 
traditional educational models, such as virtual education..., were born, have been making 
it possible for new education providers to carry out educational work more smoothly than 
ever before. The emergence of the market for higher education services in the globalization 
trend is increasing [4]. Many university programs cross borders from universities through 
the internet and other means of distance education to learners worldwide. One of the new 
trends is distance education through “virtual universities” that can form and develop. A 
form of transition from an real university to a Metaverse university is the Hybrid university 
(Augmentation model — replicates everything and enables access anytime and form 
anywhere — using IT to virtualize the real environment as a means of broadening access [5]). 
IT has become a driving force in the formation of a global academic environment, and it also 
contributes to the connection, dissemination, and exchange of all achievements between 
universities, research institutes, production facilities and corporate groups around the world. 
There are many types of universities that apply information technology in training. However, 
main form of deployment is via the internet with e-learning. Currently, some popular forms of 
training by e-learning are: Technology - Based Training; Computer - Based Training; Web- 
Based Training; Online Learning/Training; Distance Learning. Nowadays, training model in 
the Metaverse - Metaversity environment has been formed a in the world. 


DEFINITION 


Metaverse 


Neal Stephenson coined the term “metaverse” for the first time in his popular novel 
“Snow Crash” in 1992; and then Facebook rebranded as Meta and declared the metaverse 
to be Silicon Valley's “next big thing” in 2021 [1]. 

The Metaverse is the post-reality universe, a perpetual and persistent multiuser 
environment merging physical reality with digital virtuality. Itis based on the convergence of 
technologies that enable multisensory interactions with virtual environments, digital objects 
and people such as virtual reality (VR) and augmented reality (AR). Hence, the Metaverse is 
an interconnected web of social, networked immersive environments in persistent multiuser 
platforms. It enables seamless embodied user communication in real-time and dynamic 
interactions with digital artifacts. lts first iteration was a web of virtual worlds where avatars 
were able to teleport among them. The contemporary iteration of the Metaverse features 
social, immersive VR platforms compatible with massive multiplayer online video games, 


open game worlds and AR collaborative spaces [12]. 
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Metaversity 


A university inside the Metaverse [4]. Itis a concept based on combining Metaverse 
and University. It is a Metaverse community-based education platform focused on 
empowering individuals with the skills and knowledge needed to teach and learn, as well 
as activities related to education and training. In addition, metaversity also helps training 
institutions and educators... generate income on the web3 in the Metaverse environment. 
The application of Metaverse in training implementation has been forming a “market” for 
cross-border education and a sustainable creative industry. Multimedia and interactive 
display in Metaversity is a new display method that gives learners a variety of experiences, 
excitement, and fun besides receiving one-way information. They can actively interact with 
each other, collaborate, and expand to learn more information and knowledge in all relevant 
fields in the metaverse [1]. 

The Deputy Minister of Information and Communications affirmed during the 
conference “The future of the internet”, over the past 20 years, Vietnam's internet has 
developed strongly with a usage rate of 70.3%, becoming an essential need, and the 
foundation for the development of the digital economy and digital society. This is a great 
potential when digital transformation in higher education and virtual university application 
in Metaversity. If metaverse is applied in Vietnamese education, it will be a huge change, 
bringing practical benefits to learners, lecturers, and administrators, helping to access 
necessary information for managing and teaching and learning anywhere, anytime with good 
experiences in metaversity. The university will be better managed, more transparent, more 
economical and the quality of education will be improved. For examples: Greenwich Vietnam, 
also recently launched a new campus on Metavert (https://meta.horizonland.app/) [4]. It was 
only integrated into Metarvert at first for enrollment and communication purposes. However, 
this has made a worthwhile contribution to increasing its reputation with the Greenwich 
Vietnam into community and give learners a good experience in Metavert's ecosystem.The 
implementation of integrated teaching and learning on metaverse will be expanded and 
create a large educational ecosystem of Metaversity. With a virtual classroom supported 
by applications in Metaverse, participating learners can conduct learning and research at 
their own pace with mentor in the metaverse space; or even there is no mentor at all, simply 
learners participate in metaversity to experience, interact and create new content... There 
are virtual classrooms that contain a variety of ready-made learning materials that students 
can use. Students can follow up without the help of the instructor/ mentor. E-lecturer in this 
virtual classroom model replaces the lecturer to help learners based on pre-programmed 
situations and Al learned during e-learning operation. 

Learners can have fun and experience with excitement in class instead of just 
passively studying in physics class. Integrating the interactive display of course content 
into the Metaverse environment is a new way of applying technology that changes the way 
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a classroom is run. There are many ways to participate in learning on Metaverse. It can 
be combined with a traditional classroom at a training facility with the metaverse. It can be 
combined with a traditional classroom at a training facility with the metaverse. Learners have 
the feeling that they are really in the classroom with the best experience for learning and 
developing the necessary skills. It provides the ability to visualize complex phenomena and 
gather information in the clearest and most complete way. In the metaversity environment, 
the learner themselves who will decide how and with whom to study... However, to maintain 
the quality of metaversity, it is necessary to apply the same process management and 
quality management philosophies to ensure a good training system, a good trainer system 
and a good digital content system. Metaversity creates an open higher education system 
with open input, open training, open assessment from many sides, and many perspectives. 


METAVERSITY MODEL 


Metaversity is a digitally built university using VR in the metaverse; a digital university 
campus is both in sync with and out of sync with the real world. Live Classes were taught 
by a qualified lecturer. And his lecture can be recorded and deployed to the content library 
in Metaversity for students to visit and experience for themselves later. Metaversity facilities 
make physical space “flat” (geographically), meaning that learners and teachers will interact 
in a virtual classroom environment in metaversity, even though they may be far apart. Smart 
devices are installed in the classroom to digitize the physical version of the classroom and 
exist in the Metaverse space as an extension and enhancement of the learner's experience 
[3]. In the near future, Metaverse as a large society, combined with many technology 
businesses, also creates conditions for students to participate in virtual reality experiences 
and have the opportunity to access high technology. On the basis of this research, the paper 
proposes the following Metaversity model: 


Metavergse 


Ny 
S o 
7 
Edu Experience 
Metaverse 
Metaversity 


Main Metaversity elements 


Main Metaverse elements 


Metavergse 


Fig. 1. Metaversity model 
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Metaversity is the result of integrating an entire university into Metaverse: Human 
(lecturers, administrators, students...); curriculums and teaching, learning materials; 
teaching methods; assets, facilities... Allelements in the real world are digitized and operated 
in the Metaversity environment to provide the best teaching and learning experience 
(administrators, lecturers, students...) 


Learning and teaching ecosystem 


Metaversity (learning and teaching ecosystem in the metaverse) emerges and 
creates a social space for people to learn and share experiences in productive ways: 
Well-designed immersive environment with interactive way of visualization makes learning 
sessions focused and engaging with stronger and more memorable visuals. They can be 
created in the metaverse, both passive learning (the teacher conveys knowledge) and 
comprehension through reading (the learner read) is replaced by truly experiential learning. 
Learners can grab, manipulate, and walk around learning objects; they can be immersed in 
historical events; among many other opportunities. 

In education and training, the learning-by-doing process is the best way to learn 
to accumulate knowledge and form skills, even if itis a large and complex project. As a 
result, learners are “embedded” in an immersive metaverse while they also use external 
peripherals to further enhance the learning experience. Overall, the metaverse is a very 
safe place for learners to “just be themselves” so that they actively learn at their own pace, 
making mistakes and making mistakes again, until the learning goal is accomplished. 
Moreover, learning from mistakes is also a good way to master skills. If learners do this in 
the real world, it will often be very expensive or impossible to do. 


Metaverse also removes spatial barriers 


No matter where learners are, they can enter the immersive world for any learning 
purpose. Usually, they are spatially sustainable. The learners and teachers from all over the 
world can interact and work as a team in the metaverse based on a conceptual model just 
like the real world (both in space and time). 

The most compelling feature of the metaverse is that learners can interact with 
each other and with the instructor, no matter where they are located. Immersive Metaverse 
enables social and community learning experiences. Today, according to the International 
Telecommunication Union, 369 million young people do not have access to the Internet, 
resulting in exclusion, fewer educational resources and limited opportunities for children 
and young people to access the Internet. UNESCO also reports that 260 million children 
between the ages of 5 and 16 are out of school. Metaverse will be a connected ecosystem 
that brings opportunities and access to quality education for all. It acts as an educational 
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bridge in a very effective way. Universities need to work with customers, governments,... 
to bridge the digital divide in society, and promote educational programs like Connect To 
Learn, which help achieve the goals and Global Sustainable Development Goals at the 
same time. The transition from offline “classroom” to “Metaverse”, this is truly a journey to 
democratize education [1]. 

Due to the better interaction in the spatial and social structure of the metaverse, 
lecturer are also more likely to make requests to enhance work. A properly designed 
Metaversity environment will allow for the creation of entirely new ways of teaching and 
learning. Instructors can use role-playing techniques with avatars to present real-life 
situations. Learners can perform assignments easily because the 3D environment in the 
metaverse allows to simulate reality accurately [8]. From there, metaversity will build a skill 
database of learners (collecting different skills and storing them in a standardized skills 
database). 

Learners can chat interactively and lively. In addition to video lectures, hosting live 
chats can be a good way to engage the learning community and make the learning process 
a lot more dynamic. A good learning management system will provide ways to interact with 
students in a virtual classroom. Interactions can be made via video, audio systems or can 
be live chat based on text input in the interactive screen. Thereby, it is possible to connect 
everyone in the virtual classroom together to enhance the student's learning experience [9]. 

Incorporate instant feedback in Metaversity. Everyone feels encouraged to explore 
with instant feedback on how well they are doing. This can be done in virtual classrooms 
by taking advantage of one of the most popular features of a learning management system: 
automated assessment. These automated assessments will map to the e-lecturer's 
expressions and actions so that students can see the results of their activities in the virtual 
classroom. 

Metaverse helps deploy e-learning training anytime, anywhere, imparting knowledge 
on demand, providing quick feedback [9]. Learners can access the courses anywhere such 
as at the office, at home, at public Internet spots, 24 hours a day, 7 days a week. Join their 
virtual classroom on metaverse to meet the e-lecturer anytime they log-in to the virtual 
classroom. 

Metaverse makes Virtual Classrooms more accessible to learners. Learners will no 
longer have to wait until class time to ask their lecturer questions: in a virtual classroom, 
they can be easily contacted to solve any problems while they are practicing [3]. Virtual 
classrooms with ready-made databases, role-playing-like organization... 

Metaverse makes e-learning deployment more flexible. Students can choose from 
online instructor or self-interactive courses, adjust their learning pace according to their 
ability, and can enhance their knowledge through online libraries and supporting of e-lecturer 


[2]. 
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Applying Metaverse 


Applying Metaverse can focus on two main contents: in educational management 
and in teaching, learning, testing and evaluation: (1) In educational management, including 
digitizing management information, creating new system of large interconnected databases, 
application of Al technologies, blockchain, data analysis,...) to manage, operate, forecast, 
support decision making quickly and accurately; (2) In teaching, learning, testing and 
evaluation, including digital materials (e-textbooks, electronic lectures, e-learning lecture 
warehouses, multiple-choice question banks), digital libraries, laboratories,... virtual 
experience, deploy online training system, build virtual universities (cyber universities). 

To build e-learning in the Metaverse environment, it will be necessary to convert 
traditional classrooms to digitize — visual/audio digital exercises for 24/7 learning and 
interaction. Learners and their mentors will automatically rely on the digital data to design a 
training program that is suitable for their ability and level. Lectures will be based primarily on 
supportive interactions and intensive teaching based on digitized fundamentals. Metaverse 
supports content innovation, teaching and learning methods, testing, and evaluation. 
Support effective implementation of synchronous solutions (including digital data e-learning 
lectures, electronic lesson design software, simulation software and teaching software). loT, 
big data, artificial intelligence and machine learning AV/VR, security, blockchain, chatbots, 
Increased Accessibility... build and form a unified training ecosystem to connect hardware 
with software into smart, energy-efficient products and easily expand the system and 
Metaverse space. Grasping this trend of new digital technologies, universities can convert 
lectures into videos, play and learn (Gamification) to virtual reality (VR) and AR applications. 

The simplest way is to use AR to add explanatory information about the images in the 
real space of the classroom, or simply on each image section of the reference book or a part 
of the classroom space with virtual classroom space in Metaverse. The Integrated Hybrid 
Classroom (combine online and offline at the same time) can even use AR to render 3D 
(digital versions) of models. Universities will not need to build large classrooms /lecture halls 
with a capacity of thousands of students. Instead, the university just builds presentations 
that can be created and saved in a virtual environment using digital technology or can be 
uploaded to Metaverse, setting up an online classroom directly in Metaverse... Learners 
can easily enter the classroom, access the lecture content data, and can download each 
part to leam (based on compliance with copyright law). 

A lecturer's workflow in Metaverse could be: (1) first creating original immersive 
content that can be hosted in metaversity; (2) then this content will be assigned to a 
blockchain such as NFT and ensure that the blockchain used by the NFT is the same as 
used by the metaverse; (3)publish content in metaversity, assign rights and avail of viewing, 
using, modification and possible resale of the content [1]. Metaverse allows content creators 
with blockchain technology to turn into assets (cryptocurrency). These digital assets can be 
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used in the real world through rendering devices or traded under NFT as is the case with 
cryptocurrencies. 

Even learners become content creators and can leverage their “property”. A student 
might create a product as part of a project assignment. Once converted and become part 
of the blockchain, that “creative” product will allow students to monetize the creation of 
digital assets. Thereby, learners can use their metaverse to communicate, connect, play 
some games, exhibit important solutions, co-create business projects, etc. Metaverse has 
gradually improved in terms of platform, scale, security, and partnerships to achieve one 
result: a dramatically enhanced learning experience that improves learning outcomes and 


the teaching experience. 


Disadvantages of e-learning on Metaverse 


In addition to the outstanding potential of being applied in the future, e-learning on 
Metaverse also has some disadvantages. Some common problems such as: emotional 
and spatial issues that make an impression on learners; Direct interaction with learners is 
limited due to limitations of technology. Need to be connected to high-speed internet (5G) 
and connect to many network-connected devices to learn. If there is a problem with the 
equipment or the transmission network, learning will be interrupted; Haven created a lively 
space like a traditional classroom. 

Besides, Metaverse also requires computing power, and needs a large human 
resource of engineers, designers, and administrators to keep Metaverse up and working. 
Metaverse will need a large-scale technology infrastructure, from computing power to 
compute, 3D visualization to content for financial and trading systems. Going forward, 
Metaverse is still a technology built by many organizations. In addition, it is necessary 
to develop policies related to learning materials such as intellectual property, copyright; 
related to the quality of teaching in the network environment such as network information 
security [10]; related to personal information protection, information security in the network 
environment; and regulations related to the conditions for organization of online teaching 


and learning, quality accreditation, legality, and recognition of learning results... 


CONCLUSION 


In short, the digital era has now created very favorable opportunities to improve 
training quality and develop sustainably. Metaverse and gamification is one of the new 
trends in the context of the industrial revolution 4.0, helping to enhance the learner's 
experience while improving the quality of training. It helps us a lot in training management, 
in all stages of the teaching and learning process: from preparing lesson outlines, designing 
programs, lectures, implementing subjects, exchanging and interacting with people 
learning, grading and exhibition, communication activities... can all be integrated with 
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multimedia in metaverse. By integrating multimedia, technology in the teaching and learning 
process will enable rich and vivid content transmissions in real or virtual contexts. VR or 
AR can help to enhance the interpretation of complex ideas, explore human knowledge, 
increase interaction opportunities, provide user experiences... Vietnam is on the road of 
Industrialization, Modernization, deep international integration in all areas of life, including 
education and training... The research and deployment of e-learning applications in the 
metaverse environment will help Vietnamese education develop quickly and sustainably. 
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ABSTRACT: In this work, the heat transfer 
analysis of the three-phase oxygen reactor 
with a spiral baffled jacketed reactor are 
performed. The required number of oxygen 
reactors is analysed to provide enough heat 
input for different hydrogen production rates. 
Two types of fluids, which are helium gas 
and molten CuCl, are investigated to transfer 
heat from the jacket side to the process side 
of the oxygen reactor. In the analysis, the 
Cu-Cl cycle is assumed to be driven by a 
nuclear reactor where two types of nuclear 
reactors are examined as the heat source 
to the oxygen reactor. These types are 
the CANDU Super Critical Water Reactor 
(CANDU-SCWR) and High Temperature 
Gas Reactor (HTGR). In this work, it was 
found that the dominant contribution to the 
thermal resistance of the jacketed oxygen 
reactor system was from the reactor wall, 
where heat transfer occurred by conduction 
only. This contribution is about 80% of the 
total thermal resistance. It was also shown 
that a better heat transfer rate is required for 
SCWR than that for HTGR. Moreover, from 
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the study of the fluid types that can be used 
in the service side, it was recommended to 
use helium gas instead of molten CuCl as 
a heating fluid in the jacket. Finally, it was 
found that the size of the oxygen reactor 
must be specified from the heat balance 
studies rather than material balance. 
KEYWORDS:  Cu-Cl cycle; hydrogen 
production; oxygen; heat transfer; Spiral 
Baffled Jacket 


INTRODUCTION 


Hydrogen as a clean energy carrier 
is frequently identified as a major solution to 
the environmental problem of greenhouse 
gases resulting from worldwide over 
dependence on fossil fuels. It is widely 
believed that hydrogen will be a significant 
contributing factor to sustainable energy 
supply in the future, since hydrogen reduces 
pollution of the atmosphere that contributes 
to climate change by reducing greenhouse 
gas emissions. Despite the exponentially 
growing need for hydrogen, the key 
challenge facing the hydrogen economy 
is the lack of sustainable production of 


hydrogen (without dependence on fossil 
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fuels) in large capacities at lower costs than existing technologies. Nuclear energy offers 
some potential as one of the sources for large scale of sustainable production of hydrogen. 

Thermochemical cycles are the promising alternatives that could be linked with 
nuclear reactors to thermally decompose water into oxygen and hydrogen, through a 
series of intermediate reactions. Copper-chlorine (Cu-Cl) cycle is identified by the Argonne 
National Laboratories (ANL) as one of the most promising lower temperature cycles (Serban 
et. al., 2004). This cycle consists of four reactions, three thermal and one electrochemical. 
The four reaction steps of the Cu-Cl cycle are: 


4CuCl (aq) — 2CuCl, (aq) + 2Cu (s), electrolytic process at 25ºC step (1) 
2Cu(s) + 2HCI (9) — 2 CuCl (1) + H, (9), hydro-production at 450ºC step (2) 
2CuCl,(s) + H,O (9)  Cu,OCI(s) + 2HCI (9), hydrolysis at 375ºC step (3) 
Cu,OCI,(s) — 2CucCl (1) + 5 O, (9), thermolysis at 550ºC step (4) 


In the oxygen production step of the Cu-Cl cycle (step 4), an intermediate compound, 
solid copper oxychloride (Cu,OCL), is decomposed into oxygen gas and molten cuprous 
chloride (CuCl). The solid feed of anhydrous solid Cu,OCI, is supplied to the oxygen 
production reacior from the CuCIl, hydrolysis reaction (step 3) that operates at a temperature 
range of 350-450ºC. Gas species leaving the oxygen reactor include oxygen gas (which is 
evolved over a range a temperature range of 450 to 530ºC) and potentially impurities of 
products from side reactions such as CuCl vapor, chlorine gas, HCl gas (trace amount) and 
H,O vapour (trace amount). The substances exiting the reacior are molten CuCl and reactant 
particles entrained by the flow of molten CuCl. In the oxygen reactor, the decomposition of 
Cu,OCI, to oxygen and molten CuCl is an endothermic reaction requiring a reaction heat 
of 129.2 kJ/mol and a temperature of 530ºC, which is the highest temperature in the Cu- 
Cl cycle (Zamfirescu et. al., 2010). Thus, heat must be added to raise and maintain the 
temperature of the bulk inside the reactor. 

Thermochemical properties of copper oxychloride have been examined by Ikeda 
and Kaye (2008) and more recently by Trevani et al. (2011). Since copper oxychloride is 
not commercially available, methods of synthesis were developed in these two studies. 
The method adopted by Ikeda and Kaye (2008) involved the use of stoichiometric amounts 
of CuO and CuCl,. Lewis et al. (2005), described briefly the experimental status of each 
of the Cu-Cl reactions. They investigated the evolution of oxygen when the Cu,Cl,O was 
heated to 530ºC. They further calculated that the total amount of oxygen recovered was 
after 25 minutes. Zamfirescu et al. (2010) examined the relevant thermo physical properties 
of compounds of copper that are used in thermochemical water splitting cycles. They 
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identified the available experimental data for properties of copper compounds relevant to 
the Cu-Cl cycle analysis and design (Cu,OCI,, CuO, CuCl, and CuCl). They also developed 
new regression formulae to correlate the properties, which include: specific heat, enthalpy, 
entropy, Gibbs free energy, density, formation enthalpy and free energy. The properties 
were evaluated at 1 bar and a range of temperatures from ambient to 675-1000K, which 
are consistent with the operating conditions of the cycle. Marin (2012) provided new 
experimental and theoretical reference for the scale-up of a CuO*CuCI, decomposition 
reactor with consideration of the impact on the yield of the thermochemical copper-chlorine 
cycle for the generation of hydrogen. 

Abdulrahman et al. (2013a, 2016a, 2022a) have examined the scale-up feasibility 
of the oxygen reactor from the perspectives of the optimum size and number of oxygen 
reactors for different oxygen and hydrogen production rates. They specifically analysed 
the factors contributing to the oxygen reactor size. It was shown that the reactor size is 
significantly influenced by residence times, hydrogen production rate, mass and heat 
transfer. Abdulrahman (2016b, 2013b, 2019a, 2019b, 2022b) has investigated the scaleup 
analyses of the oxygen reactor from the perspective of indirect heat transfer using a half 
pipe and spiral baffled jacketed reactor and a helical tube inside the reactor (Abdulrahman, 
2016c). He has concluded that the size of the reactor calculated from the perspective of heat 
balances is more than that calculated from the perspective of material balances. Different 
experiments have been performed to examine the hydrodynamics and direct contact heat 
transfer in the oxygen reactor (Abdulrahman, 2015, 2016d, 2016e, 2018, 2020a, 20190). 
In the experiments, empirical equations have been formulated for the gas holdup and the 
direct contact heat transfer coefficient in the oxygen reactor. CFD simulations have been 
created to examine the hydrodynamics and heat transfer for the oxygen slurry bubble 
column reactor using two dimensional Eulerian-Eulerian approach (Abdulrahman, 201 6f, 
20169, 2020b, 2020c, 2022c, 2022d, 2023a, 2023b, 2023c; Nassar, 20283). 

In this work, the heat balance of the oxygen reactor in the Cu-Cl cycle is investigated 
analytically for different hydrogen production rates using a spiral baffled jacketed reactor. 
The thermal resistance of each section of the oxygen reactor system is investigated to 
examine the effect of each section on the heat balance. The heat balance of the oxygen 
reactor system is studied for Continuous Stirred Tank Reactor (CSTR) type that is heated by 
using a spiral baffled jacket. Two types of heat sources are studied which are Super Critical 
Water Reactor (SCWR) and High Temperature Gas Reactor (HTGR). Moreover, the type of 
the working fluid in the service side of the oxygen reactor system is examined for two types 
of fluids which are helium gas and CuCI molten salt, and a comparison between both fluids 
is performed. 
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THERMAL RESISTANCES 


The continuous stirred-tank reactor (CSTR) is usually used for multiphase reactions 
that have fairly high reaction rates. Reactant streams are continuously fed into the vessel, 
and product streams are withdrawn. In CSTR, heating is achieved by a number of different 
mechanisms. The most common one involves the use of a jacket surrounding the vessel. 

In a realistic continuous situation, where the reactor contents are at constant 
temperature, but with different service side inlet and outlet temperatures, the heat flow 
equation can be expressed by: 


Q=UA,ATm (1) 


where 4T,, is the log mean temperature difference between the bulk temperature 
of the reactor contenis (cold temperature), T,, and the temperature in the service side, T,, 
Since the inside reactor temperature is assumed constant (7,=530 ºC), there is no effect 
of heat transfer configuration (parallel or counter flow) in the equation of 4T,, and can be 


written as: 


The methodology of calculations depends on the comparison of thermal resistances 
of each section in the oxygen reactor with the calculated total thermal resistance of the 
reactor system (R,). For a hydrogen production rate of 1 kg/day, the total amount of heat 
required in the oxygen reactor can be calculated as follows; 


530 
Q=AH,é+n CPcu,0CI dT =870W (3) 
375 

where 4H = 129.162 kJ/mol, Fe 0.5 kmol/day, ts = 0.5 kmol/day, CPouoct, = 
134 J/mol.K (Zamfirescu, 2010). The limits of the integral are from 375ºC, which is the 
temperature of the fed solid, to 530ºC, which is the temperature of the decomposition 
process. For a hydrogen production rate of 100 ton/day, is 50000 kmol/day and Q 
will be 87 MW. 


For a high temperature gas reactor (HTGR), the nuclear reactor exit temperature 


Cup0Clp 


is about 1000ºC. The inlet temperature of the heating fluid in the jacket (same as the exit 
temperature of the intermediate heat exchanger (IHE)) is about 900ºC (Natesan, 2006). 
The exit temperature from the jacket is assumed to be 540ºC (because the decomposition 
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temperature is 530ºC). For a hydrogen production rate of 100 ton/day and number of 
reactors N, the total thermal resistance of the HTGR can be calculated as follows; 


1 AT, 
Rurer = TA, E =1.15x108N KW (4) 

For CANDU supercritical water reactor (CANDU-SCWR), where the nuclear reactor exit 
temperature is about 625ºC (Chow & Khartabil , 2007; Duffey & Leung, 2010) and the inlet and 
outlet jacket temperatures are about 600ºC and 540ºC respectively, the total thermal resistance 
is R = 3.55 x 107 NKAN for a hydrogen production rate of 100 ton/day. 


SCWwR 


HEAT TRANSFER ANALYSIS 


The resistance to the heat transfer is a composite of the resistances through the 
various sections indicated in Fig. 1. Using the classical film theory and heat conduction 
through composite layers, the total thermal resistance (R,) can be expressed as; 

R=R,+R,= Ry= Ro= Ros 


Ds 
E E e O (5) 
“hprDpHa nDpHa 2mkwHp hsmDsHa mDshHp 
The largest thermal resistance in Eq. (5) dominates the value of the overall heat 


transfer coefficient. 


Reac 


Boundar Bounda 
Proce Seni 


F 


Fig. 1 Thermal resistances through the reactor wall sections. 
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Heat transfer through reactor wall 


The stainless steels are the most frequently used corrosion resistant materials in 
the chemical industry. Types 321, 347 and 348 stainless steels are advantageous for high 
temperature service (such as in the oxygen reactor) because of their good mechanical 
properties. They continue to be employed for prolonged service in the 427ºC to 816ºC 
temperature range. The physical properties of Types 321, 347 and 348 are quite similar and, 
for all practical purposes, may be considered to be the same. 

The ASME Code provides formulas that relate the wall thickness to the diameter, 
pressure, allowable stress, and weld efficiency. Since they are theoretically sound only for 
relatively thin shells, some restrictions are placed on their application. For cylindrical shell 
under internal pressure, the thickness in inches units is; 


PRp 


"p= quo "* 


conditions: t<0.25DR, P=<0.385 SE (6) 

where tis the corrosion allowance in inches, AR, is the outside radius of the cylindrical 
shell in inches and Sis the maximum allowable working stress in psi. Pis the total pressure 
which is the sum of the static pressure and operating pressure of the oxygen reactor. The 
value of the joint efficiency, E, is between 0.6 and 1 (Abdulrahman, 2019b). In order to allow 
for possible surges in operation, it is customary to raise the maximum operating pressure 
by 10% or 0.69-1.7 bar over the maximum operation pressure, whichever is greater. The 
maximum operating pressure in turn may be taken as 1.7 bar greater than the normal 
(Couper et. al., 2005). 


Heat transfer through spiral baffled jacket 


The cross section of a spiral baffled jacket is shown in Fig. 2. The spacing between 
the jacket and reactor wall depends on the size of the reactor, however, it ranges from 
50 mm for small reactors to 300 mm for larger reactors. In heat transfer applications, this 
jacket is considered a special spiral baffling case of a helical coil if certain factors are used 
for calculating outside film coefficients. The leakage around spiral baffles is considerable, 
amounting to 35-50% of the total mass flow rate. Due to the extensive leakage, the effective 
mass flow rate mm” in the jacket is usually taken as 60% of the actual flow rate Mm to get a 
conservative film coefficient (Coker, 2001). 

abaigmedo-—++ (7) 
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Fig. 2 Cross section of spiral baffled jacket. 


where T, and T,, are the jacket inlet and outlet temperature respectively. The velocity 
of the flow in the jacket can be calculated from the effective mass flow rate by; 


m' 


GA (8) 


ur 


where A is the spiral baffled jacket cross sectional area and p is the density of the 
service side fluid (heating fluid). The pitch of the spiral baffled jacket can be calculated from; 


a (9) 
P=y 


where w is the width of the spiral baffled jacket. The number of coils of the spiral 
baffled jacket can be calculated from; 


n=— (10) 


where H, is the height of the reactor. The equivalent heat transfer diameter for the 
spiral baffled jacket is (Coker, 2001); 
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At a given Reynolds number, heat transfer coefficients of coils, particularly with 
turbulent flow, are higher than those of long, straight pipes, due to friction. This also applies 
to flow through a spiral baffled jacket. The equation of heat transfer should be multiplied by 
a turbulent flow correction factor, involving the equivalent diameter and the diameter of the 
spiral coil. This correction factor is equal to (1 + 3.5 (D /D)), where D, is the equivalent heat 
transfer diameter for the spiral baffled jacket, which is calculated from (Coker, 2001): D = 4 
w, and Dis defined as the centerline diameter of the jacket passage. 

At Re > 10,000 the Sieder-Tate equation for straight pipe can be used to calculate the 
outside film coefficient (Coker, 2001); 


0.14 D 
Nu = 0.027 Reº8 Prº33 (=) f1 435 o) (12) 
Uw Dc 
where Dis defined as the centerline diameter of the jacket passage and is calculated 
as; 


Do —- D: 
De = Dj + DA 


(13) 
where D, and Do are the inlet and outlet jacket diameter respectively. At Re < 2,100, 
the following equation can be used (Coker, 2001); 


0.14 ,D 0:33 
Nu = 1.86 Reº33 pr033 (02) (5) (14) 
Ww 


where L is the length of the spiral baffled jacket. For 2,100 < Re < 10,000, above 
equations can be used depending on the value of Re (Coker, 2001). 


Process side heat transfer 


Oxygen reactor is a multiphase reactor that contains solid particles (Cu,OCI,), 
molten salt (CuCl) and oxygen gas (O,). In the study of indirect heat transfer, the presence 
of oxygen gas is neglected, because it is assumed that the oxygen gas will leave the reactor 
immediately after it is formed. The presence of solid particles will affect the viscosity of the 
molten salt. Since the density of the Cu,OCI, solid (4080 kg/m?) has a value near to that 
of the CuCl molten salt (3692 kg/mº), it is assumed that solid particles and molten salt are 
well mixed and form homogeneous slurry. This well mixed homogeneous slurry will lead to 
a more uniform temperature profile inside the oxygen reactor, that's why the temperature 
profile is assumed to be constant and equal to 530ºC. Thermo physical properties are 
calculated for the slurry mixture at 530ºC. 
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The dynamic viscosity of slurry can be described as relative to the viscosity of the 
liquid phase; 


Ugn=H MH (15) 


where u, is the relative dynamic viscosity (dimensionless), u, is the dynamic viscosity 
of the liquid (CuClI molten salt). Depending on the size and concentration of the solid 
particles, several models exist to describe the relative viscosity as a function of volume 
fraction 6 of solid particles. 


V. 
g=—— 


= (16) 
Vszr 


where Vç and V,,, are the volumes of solid particles and slurry respectively. In the 


case of extremely low concentrations of fine particles, Einstein's equation (Einstein, 1906) 
may be used; 


n=1+25 6 (17) 


In the case of higher concentrations, a modified equation was proposed by Guth and 
Simba (1936), which considers interaction between the solid particles; 


H,=1+25 0+141 6º (18) 
Other thermo physical properties of the slurry mixture can be calculated from the 


volume percent of the solid and molten salt. For example, the average density of the slurry 
can be calculated as follows; 


Pan= Ps s+ Pp (1-O9 (19) 


where pç and p, are the densities of solid and liquid respectively. The average specific 
heat of the slurry is; 


C + p, Cp, (1 — 
Chai Ps Cps bs + p, Cp ( Ds) (20) 
PsLR 
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where Cp, and Cp, are the specific heats of the solid and liquid respectively. A first 
order estimate of the effective thermal conductivity of a fluid filled porous media can be 
made by simply accounting for the volume fraction of each substance, giving the resulting 
relation based on the porosity and the thermal conductivity of each substance (Guth & 
Simba, 1936). Hence, the effective thermal conductivity of the slurry is calculated from; 


kan=k:05+k (1-05) (21) 


where kç and k, are the thermal conductivities of the solid and liquid respectively. For 
an agitated reactor, the inside film heat transfer coefficient (A,) can be calculated from the 
following Nusselt number correlation (Coker, 2001); 


c 
Nup = CReº Pr! (57) (22) 
w 


where, Cis constant, Re = Reynolds number = — and the agitator diameter D, 
= D,/3. The values of constant C and the indices a, b and c depend on the type of agitator, 
the use of baffles, and whether the transfer is to the vessel wall or to coils. An agitator is 
selected on the basis of material properties and the processing required. The heat transfer 
forms part of a process operation such as suspending or decomposing solid particles and 
dispersing the oxygen gas in the molten salt. Several methods for selecting an agitator are 
available (Kline, 1965; Himmelblau et. al., 2008). Penny (1970) showed one method based 
on liquid viscosity and vessel volume. According to the Penny's graph and the specifications 
of oxygen reactor contents and size, the type of the impeller that can be used is a propeller 
with 420 rpm. For baffled reactor, three blades propeller and transfer to reactor wall, the 


Nusselt number equation is (Coker, 2001); 
(Re > 5000) (23) 


Up po 


Nup = 0.64Reº$” prº.33 ( 
Uw 


Fouling resistance 


An important part of the specification for the oxygen reactor is the assignment of 
fouling effects. A recommended way to provide an allowance for fouling, is the use of 
individual fouling resistance values R, and R,, for the two sides of the oxygen reactor 
as used in Eq. (5). R, is the fouling resistance that occurs on the internal surface of the 
reactor (process side) as a result of deposits that accumulate from the reactor contents. 
Ro is the fouling resistance on the service side of the reactor (like jacket). For helium gas 
heating fluid, R,, is expected to be negligible since there should be no build-up associated 
with clean dry helium. Values of the fouling resistances are specified which are intended to 
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reflect the values at the point in time just before the exchanger is to be cleaned. Fouling of 
heat exchange surfaces can cause a dramatic reduction in the performance of the reactor 
because of the relatively low thermal conductivity of the fouling material. 

There are different approaches to provide an allowance for anticipated fouling in the 
design of oxygen reactor. In all, the result is to provide added heat transfer surface. This 
generally means that the reactor is oversized for clean operation and barely adequate for 
conditions just before it should be cleaned. 


TYPE OF WORKING FLUID IN THE SERVICE SIDE OF OXYGEN REACTOR 


Two types of fluids are highly recommended as working fluid in the oxygen reactor; 
CuCI molten salt and high-pressure helium gas (He). Helium gas is recommended amongst 
the other noble gases because of its cnemical inertness and the relatively good transport 
properties (Rousseau & Van Ravenswaay, 2003; Wang et. al., 2002; Kikstra & Verkooijen, 
2000). 

There are substantial differences between molten salts and high-pressure helium 
that must be considered in selecting the working fluid as a heating medium in the oxygen 
reactor. These key differences are thermal performance, materials compatibility, and safety. 


Thermal Performance 


The thermo-physical properties for high-pressure helium and CuCI molten salt are 
summarized in Table 1. 


Material scan SM) Helium (7.5 MPa) Helium (2 MPa) 

T menina (00) 430 ; : 

T posa (0) 1490 ; : 

p (kg/m? 3692 35 0.939 

C, (kJ/kg.ºC) 0.66 52 5.2 

PC, (kJ/mº.2C) 2450 182 4.883 

k (W/m.ºC) 0.23 0.37 0.3687 

4 x 10º (pa.s) 260 4.7 4.7 

vx 108º (m?/s) 0.7 13.4 50 

Pr 4.29 0.66 0.6624 


Table 1 Comparison of thermo-physical properties of helium and CuCl molten salt at an average 
temperature of 750ºC 
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From Table 1, it can be seen that the volumetric heat capacity, p C,, of molten 
salt is over two orders of magnitude greater than that of high-pressure helium. The much 
higher p Cof molten salts, compared to high-pressure helium, has a significant effect upon 
the relative heat transfer capability. In general, a molten salt loop uses piping of smaller 
diameters, and less pumping power than those required for high-pressure helium. These 
differences in pumping power and pipe size reduce the capital cost of the piping system, 
and allow the arrangement of process equipment to be optimized more easily since process 
heat can be delivered over larger distances easily. Also, heat transfer coefficients for molten 
salts are typically greater than those for helium. 


Materials 


Two primary aspects must be considered in selecting the high-temperature materials 
and heating fluid used for the oxygen reactor: corrosion and high-temperature mechanical 
properties (strength, creep, and fabric ability). Clean helium clearly does not have the 
potential to corrode loop materials, while, molten salts exhibit higher corrosion rates. The 
potential thermal performance advantages of molten salts suggest that the high-temperature 
corrosion testing with molten salts should be a priority for the oxygen reactor design. 


Safety 


The copper-chlorine process uses chemicals which are hazardous. The heat transfer 
fluid in the service side of oxygen reactor provides a stored energy source that can potentially 
be released rapidly, generating mechanical damage and potentially dispersing flammable 
or toxic chemicals. For helium gas, the stored energy comes from the high pressure of 
the gas, and the large volume of gas due to the large duct sizes required for transferring 
helium with reasonable pressure losses. For molten salts, the stored energy comes from 
the high temperature and high heat capacity of the liquid. This energy can be released if 
the molten salt mixes with a volatile liquid (e.g. water), through a well-studied phenomenon 
typically referred to as a “steam explosion.” Molten salts, due to their much lower pumping 
power and small piping size, permit greater physical separation between the nuclear reactor 
and the hydrogen production plant. This reduces or eliminates the need for berms or other 
structures to provide isolation between the reactor and hydrogen plant. 
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RESULTS AND DISCUSSION 


Heat transfer calculations of the oxygen reactor system 

The physical properties of the heating fluid in the service side of oxygen reactor are 
calculated at the mean temperature of the fluid, which is for a HTGR is 720ºC and for SCWR 
is 570ºC. Table 2 shows the physical properties of both helium and molten salt CuCl for both 
HTGR and SCWR. 


u 
i k w 
cas Fluid | P(MPa) | »p e ( kg ) EG Pr | Reference 
m.s 
ATOR i (Petersen, 
(720ºC) Helium 1 0.48 4.6 0.36 0.66 1970) 
SCWR ; (Petersen, 
(570ºC) Helium 1 0.48 4.13 0.32 0.66 1970) 
HTGR (Zamfirescu, 
(720º0) CuCl 0.1 3692 146 0.19 5 2010) 
SCWR (Zamfirescu, 
(570ºC) CuCl 0.1 3692 188 0.2 6.11 2010) 


Table 2 Physical properties of helium gas and CuCI molten salt for different mean temperatures 


For a helium gas with a pressure of 1 MPa and a temperature range of 570-720ºC, 
the compressibility factor range is 1.001426-1.009594 (Petersen, 1970). This factor varies 
from unity by less than 1% which means that helium gas can be regarded as an ideal gas. 
The equation of sound for gases is; 


c=vVyRT (24) 


where y is an adiabatic index (y,, = 1.6667) and Ris the gas constant (R,, = 2077 J/ 
Kg. K). By using Ea. (24), the theoretical maximum speed of a helium gas (which is 1/3 of 
speed of sound) for a HTGR (T=720ºC) is 618 m/s and for a SCWR (T=570ºC) is equal to 
570 m/s. In this work, to be more conservative for getting incompressible flow, the operating 
speed of the helium gas is considered to be 300 m/s in baffled jacket reactor for both HTGR 
and SCWR. The operating speed of the CuCl molten salt is assumed to be 3 m/s for baffled 
jacket. These values are within the range of molten salt nuclear reactors reported (Huntley 
& Silverman, 1976; Caire & Roure, 2007). 

The physical properties of the slurry mixture inside the oxygen reactor (process side) 
are calculated from Eqg. (15) to (21). Table 3 shows the physical properties of both the 
solid particles and the molten salt inside the oxygen reactor at a temperature of 530ºC 
(Zamfirescu, 2010; Marin, 2012). 
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J W Kg 

Cp Kg.K x m.K Pa 

Solid Cu,OCI, 623.7 0.451 4080 
Molten salt CuCl 650.8 0.2 3692 


Table 3 Physical properties of Cu,OCI, solid and CuCI molten salt at a temperature of 530ºC 


Heat transfer by jacketed reactor 


To calculate the thermal resistance of the reactor wall, it is necessary to know the 
thickness of this wall. The thickness of the reactor wall (t) is calculated from Eg. (6). The 
pressure, Pin Eq. (6) is the design pressure of the oxygen reactor which is taken here as 
1.7 bar greater than the normal total pressure. The normal total pressure is the sum of the 
operating pressure (P,) that is equal to 1 bar (14.5 psi) and the static pressure (P9). Thus 
the design pressure is; 


P=P+P+1.7(ba)=pgH,+P+1.7 (bar (25) 


where p is the density of the slurry that is calculated from Eqg. (19). The size of 
oxygen reactor has been estimated on the basis of the hydrogen production scale, mass 
balance, residence time, and aspect ratio, among others. The diameter range of 3-4 m and 
the aspect ratio of 2 were recommended for an industrial hydrogen production scale of 
over 100 tons/day (Abdulranman, 2013a). In this work, the value of the reactor diameter is 
assumed to be 4 m and the height is 8 m. 

The material ofthe reactor wall is assumed to be stainless steel 321 and the maximum 
allowable working stress (S) for this stainless steel is 3600 psi at T=649ºC (Rozic, 1973). 
The value of the thermal conductivity (K,) for stainless steel 321 in the temperature range of 
20-500ºC is 21.4 W/m.K (Wang et. al., 2002). For known corrosive conditions, the corrosion 
allowance (t,) is 8.9 mm (0.35 in.) (Chenoweih, 1990). This value is used for the oxygen 
reactor because of the relative high corrosion susceptibility of CuCI molten salt. The joint 
efficiency E is taken as 0.8. By substituting the above values into Eq. (6), the thickness of 
the reactor with a diameter of 4 m and a heightof8 mwillbe, t=2.7 in=7 cm 

In order to calculate the thermal resistance due to fouling in the oxygen reactor, the 
fouling factor must be known or estimated. According to Tubular Exchanger Manufacturers 
Association (TEMA), the fouling factor for molten heat transfer salts is equal to 0.000088 m? 
K/W (Chenoweth, 1990). In this work, to be more conservative, twice of this value, 0.000176 
m2 K/W, is used to be the fouling factor for molten salt CuCl. 

To calculate the thermal resistance of the jacket side, it is assumed that the thickness 
of the jacket for spiral baffled jacketis D,/15, where D, is the reactor diameter. Table 4 shows 
the values of the numbers of oxygen reactors required for each section of the oxygen reactor 
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for different heating fluids and different nuclear reactors. Figure 3 shows a comparison of 
the number of reactors required between HTGR and SCWR. This comparison is for each 
section of the oxygen reactor system with a helium gas at a pressure of MPa. From Table 
4 and Fig. 3, itcan be seen that the maximum number of reactors required comes from the 
wall of the reactor. This is because of the large thickness and small conductivity of the wall 
material. The large diameter requires thick wall to provide enough mechanical stresses and 
the large height produces high static pressure inside the reactor, which in turn requires more 
thickness for the wall. It can be seen also that the number of oxygen reactors for SCWR 
are higher than that for HTGR by more than three times. This is because of the higher 
temperature difference between the service and process sides in the HTGR. This means 
that using HTGR is more efficient for producing heat to the oxygen reactor than SCWR. 
Figure 4 shows the total number of jacketed oxygen reactors with each nuclear reactor for 
both helium gas and CuCI molten salt. From this figure, it can be seen that there is no big 
differences between helium gas and molten CuCl in heating the jacketed oxygen reactor. 

Fig. 5 shows the comparison of number of oxygen reactors versus the hydrogen 
production rate between material balance and heat balances for both HTGR and SCWR 
for a residence time of 0.5 hr and for helium gas service fluid. From Fig. 5, it can be seen 
that the numbers of reactors calculated from heat balance are 18 times higher than that 
calculated from material balance for HTGR and 57 times for SCWR. That means the design 
of oxygen reactor size is controlled mainly by the heat balance. 


Nuclear Fluid 
Reactor Type Type P (MPa) N, Neo Ny Nes Novato Notarbatrte 
He 1 0 3.4 35 
HTGR 1.4 1.5 27.8 
CuCl 0.1 1.5 1.4 34 
He 1 0 13 113 
SCWR 4.9 4.9 90 
CuCl 0.1 4.8 4.6 110 


Table 4 Number of reactors of each section of the jacketed oxygen reactor system for different heating 
fluids and different nuclear reactors 
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Fig. 3 Number of reactors for each section in oxygen reactor system heated by 1Mpa helium gas for 


both HTGR and SCWR. 
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Fig. 4 Total number of oxygen reactors with each nuclear reactor for both 1MPa helium gas and CuCl 


molten salt. 
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Fig. 5 Number of oxygen reactors calculated by material balance and heat balances of both HTGR and 


SCWR versus hydrogen production rate for a residence time of 0.5 hr. 
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CONCLUSIONS 


In this work, it can be concluded that the dominant heat transfer thermal resistance 
of the jacketed oxygen reactor system, is from the reactor wall which contributes in about 
80% of the total thermal resistance for both types of heat sources (SCWR and HTGR). The 
service side contributes in about 10% and the process side and fouling in about 4% for 
each. The high thermal resistance of the wall is due to the large thickness and small thermal 
conductivity of the reactor wall. 

Heat transfer calculations showed that the total thermal resistance needed for SCWR 
is less than that for HTGR. That means a better heat transfer rate is required for SCWR 
than that for HTGR. Thus, HTGR is more efficient than SCWR in providing the required heat 
for the thermal decomposition process inside the oxygen reactor, because of the highest 
temperature difference between the service and process side in HTGR than in SCWR. 

From the comparison study of the type of fluid in the service side, it is concluded that 
the thermal performance of the molten CuCl is better than that of helium gas which in turn is 
better than molten CuCl from the perspectives of material and safety aspects. By comparing 
the number of reactors for both fluids, it is shown that there is no significant difference 
between both fluids and in this case, itis recommended to use helium gas as a heating fluid 
in the service side of the oxygen reactor system. 

From the comparison of the number of reactors obtained from material and heat 
balances, it is concluded that the size of the oxygen reactor is specified from the heat 
balance study rather than material balance. 
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ABSTRACT: The thermochemical water 
splitting process using the copper-chlorine 
(Cu-Cl) cycle is an innovative method for 
generating hydrogen gas. This process 
involves a thermolysis reaction where 
molten salt, a solid reactant, and gaseous 
oxygen interact to produce oxygen gas. 
The primary aim of this work is to contribute 
to the development and enhancement 
of the thermochemical Cu-Cl cycle for 
hydrogen production. To achieve this goal, 
a comprehensive review of the multiphase 
oxygen generation reactor, which is a crucial 
component of the Cu-Cl cycle, is conducted. 
The review covers various aspects of 
the oxygen reactor, including; type and 
description of the Oxygen reactor, material 
features, variables affecting reactor's 
size, and scale-up evaluations. This work 
aims to provide a thorough understanding 
of the oxygen reactor in the Cu-Cl 
cycle, highlighting its significance in the 
thermochemical production of hydrogen and 
the key considerations for its optimization 
and scale-up. 
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INTRODUCTION 


Hydrogen is increasingly recognized 
as a pivotal solution to the environmental 
challenges posed by greenhouse gas 
emissions, which are largely attributed 
to the global dependence on fossil fuels. 
kt is widely regarded as a significant 
contributor to the future sustainable energy 
supply, as its utilization can substantially 
reduce atmospheric pollution and mitigate 
climate change by lowering greenhouse 
gas emissions (Forsberg, 2007). Given 
its versatile applications across various 
industries, the demand for hydrogen 
is projected to rise considerably in the 
forthcoming decade. 

The majority of hydrogen production 
on a global scale currently relies on 
methods such as steam-methane reforming 
or partial oxidation of heavy hydrocarbons. 
While these processes are efficient, they 


generate substantial quantities of carbon 
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dioxide, contributing to greenhouse gas emissions. These conventional methods are facing 
increasing challenges due to escalating greenhouse gas emissions, reducing natural gas 
resources, and rising costs associated with carbon capture and storage. 

In light of the exponential growth in demand for hydrogen, the primary challenge for 
the hydrogen economy is the sustainable generation of hydrogen on a large scale, without 
dependence on fossil fuels, and at costs lower than those of existing technologies. In this 
context, nuclear energy emerges as a potential source for large-scale, sustainable hydrogen 
production. 

Thermochemical cycles represent a promising alternative that can be integrated with 
nuclear reactors to thermally decompose water into oxygen and hydrogen through a series 
of intermediate processes. These cycles facilitate the transfer of heat between different 
reactors, some of which are endothermic and others exothermic, via heat exchangers 
that either provide or recover heat from specific processes. The copper-chlorine (Cu-Cl) 
cycle, identified by Argonne National Laboratories (ANL) as one of the most promising low- 
temperature cycles, stands out in this regard (Lewis et. al 2003; Serban et. al., 2004). One 
of the most notable advantages of the Cu-Cl cycle is its lower temperature requirement 
compared to many other cycles, which enhances its efficiency and accessibility. Furthermore, 
the cycle can be conducted in a manner that necessitates minimal high-quality energy and 
solid exchanges, making it more environmentally and economically viable. 

The integration of this cycle with nuclear energy offers a sustainable pathway for 
hydrogen generation, addressing the pressing need for clean energy solutions. As the world 
moves towards a more sustainable energy future, the development and optimization of 
thermochemical cycles like the Cu-Cl cycle, in conjunction with advancements in nuclear 
technology, hold great promise for meeting the growing demand for hydrogen while mitigating 
the environmental impact of energy production. This cycle comprises three reactions: two 
thermal and one electrochemical. The Cu-Cl cycle comprises three reactions, two thermal 
and one electrochemical: 


CuCl (a) + 2HCI (9) > CuCl, (a) + H, (9) 100ºC step (1) (1) 
2Cucll, (s) + H,O (9) « Cu,OCI, (s) + 2HCI (9) 375ºC step (2) (2) 
Cu,OCI, (s) > 2CuCl (|) + 42 O, (9) 5B30ºC step (3) (3) 


where a, s, /and g denote to aqueous, solid, liquid and gas respectively. 

In the oxygen generation stage (step 3), an intermediate chemical known as 
solid copper oxychloride (Cu,OCI,) undergoes decomposition, resulting in the formation 
of oxygen gas and molten cuprous chloride (CuCl). This stage is crucial for the overall 
efficiency of the Cu-Cl cycle. The oxygen generation reactor receives an input of anhydrous 
solid Cu,OCI, from the preceding step in the cycle, which is the hydrolysis of CuCI2. This 
hydrolysis process operates at temperatures ranging from 350 to 450ºC. Once inside the 
oxygen reacior, the solid Cu,OCI, is subjected to higher temperatures, typically between 
450 and 530ºC, which triggers its decomposition into oxygen gas and molten CuCl. 
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However, the process is not without its challenges. The gaseous output from the 
oxygen reactor may include not only the desired oxygen gas but also potential contaminants 
from side reactions. These could include vaporized CuCl, chlorine gas, and trace amounts of 
hydrochiloric acid (HCl) gas and water vapor (H,0). These by-products need to be carefully 
managed to ensure the purity of the produced oxygen and the overall efficiency of the cycle. 

Moreover, the substances exiting the reactor are primarily molten CuCl, but there 
might also be solid CuCI, present. This can occur due to the incomplete breakdown of CuCl, 
at temperatures below 750ºC (Micco et. al., 2007). Additionally, reactant particles may be 
entrained in the flow of molten CuClI, further complicating the process. 

Overall, the oxygen generation stage of the Cu-Cl cycle is a complex process that 
requires careful control and management to ensure its efficiency and effectiveness in 


contributing to the sustainable production of hydrogen. 


OXYGEN REACTOR TYPE 


Three phases comprise the oxygen reactor: the solid phase (copper oxychloride 
particles), the liquid phase (molten salt), and the leaving gaseous phase (oxygen). 
In chemical, petroleum, and biological processes, multiphase reactions such as gas/ 
liquid/solid are often employed. The industry offers various different kinds of multiphase 
reactors, which may be grouped into two broad categories: fixed bed reactors and slurry 
phase reactors. Fixed bed or packed bed reactors have a stationary solid phase. Because 
the oxygen reactor's solid phase is continuous, it cannot be termed a fixed bed reactor 
(Abdulrahman, 2016a). Stirred tank reactors (STR) and slurry bubble column reactors are 
the most common continuous slurry phase reactors (SBCRs). The continuous stirred tank 
reactor (CSTR), which is often used for liquid-phase or multiphase processes with moderate 
reaction rates, is a mechanically agitated reactor in which small particles are suspended in 
the liquid phase by agitation. Continuous reaction streams are supplied into the vessel, while 
product streams are removed. In general, a continuous flow stirred tank reactor (CFSTR) 
implies that the fluid is properly mixed, ensuring that the reaction mixture's parameters 
(e.g., concentration, density, and temperature) remain uniform throughout the system. As a 
result, the conditions throughout the reactor are same, and the temperature at the reactor 
output is identical. The solid phase in the SBCR is composed of small particles suspended 
in the liquid phase as a result of gaseous bubbles supplied from the bottom, often through 
a sparger (Abdulrahman, 2016a; Li et. al., 2016). 

Oxygen reactors may be classified as CFSTRs or SBCRs. The decision between 
these two kinds of reactors is determined by the heat transfer rate efficiency, the agitation 
efficiency, the ability to scale up, and other design requirements. The oxygen reactor's 
thermal design demands adequate agitation inside the reactor and a sufficiently high rate 
of heat transfer to the solid particles. Agitation may be accomplished mechanically, as 
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in the CFSTR, or by gaseous bubbles, as in the SBCR. Although the rate of agitation is 
faster with a mechanical agitator than with bubbles, an agitator is undesirable in a highly 
corrosive medium such as an oxygen reactor due to the numerous internals required, 
such as the propeller, shaft, baffles, and other accessories that support the mixing system 
(Abdulrahman, 2016a). 

Heating may be accomplished in two ways: indirectly or directly. The two most often 
used indirect procedures are to surround the vessel with a jacket or to use an inside coil. 
Direct heating may be accomplished by the use of gaseous bubbles and is more efficient 
than indirect heating, but is more complicated to scale up. The oxygen gas that exits the 
oxygen reactor at around 530ºC may be heated to a higher temperature (such as 600ºC) 
and then re-injected into the oxygen reactor from the bottom through a sparger. Direct 
contact heat transfer between the hot oxygen gas and the slurry within the oxygen reactor 
is possible in this instance (Abdulranman, 2016a). 

According to the comparison above, SBCR is more efficient in heat transmission 
than CFSTR but more difficult to scale up. However, with SBCR, thermal energy may be 
delivered to the reactor directly through gas bubbles, which can also be utilized to agitate 
the reactor's contents. 


OXYGEN REACTOR SYSTEM DESCRIPTION 


As with many chemical reactors, the oxygen reactor needs careful regulation of heat 
transfer to operate at peak efficiency. Cu,OCI, decomposition into oxygen and molten CuCl 
is an endothermic process in the oxygen reactor, requiring a reaction heat of 129.2 kJ/mol 
and a temperature of 530ºC, the highest temperature in the Cu-Cl cycle. As a result, heat 
must be provided to increase and maintain the temperature of the reactor's bulk (Ryskamp, 
2004). The total amount of heat needed is the sum of reaction heat and the heat necessary 
to elevate the reactant temperature from 375ºC (the temperature of solid particles formed 
during the hydrolysis process) to 530ºC. The heat transfer coefficient, the difference between 
the bulk and surrounding temperatures of the heat transfer fluid (i.e., the driving force), and 
the size of the contact region where the heat transfer occurs all affect the heat transfer. 
Additionally, the heat transfer coefficient is dependent on the oxygen generation process's 
operating conditions, the physical qualities of the inventory materials, and the geometry of 
the equipment. Typically, this coefficient is established by the use of empirical correlations. 

The oxygen reaction is a high temperature reaction that requires a source of high 
temperature heat. This heat may be generated by generation IV nuclear reactor designs 
such as the sodium-cooled reactor, the Canada Deuterium Uranium (CANDU) supercritical 
water reactor (CANDU-SCWR), Canada's Generation IV reactor, and the high temperature 
gas reactor (HTGR). While the first two reactors operate at lower temperatures (510ºC 
— 625ºC), they are ideally suited for matching a low temperature Cu-Cl thermochemical 
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cycle. However, coupling low temperature cycles to the HTGR (1000ºC) may enable more 
cogeneration, resulting in much better hydrogen production efficiency. Solar energy is 
another non-polluting source of high-temperature heat. 

In an indirect cycle system, an intermediate heat exchanger (IHX) is utilized to 
transmit heat from the main fluid in the nuclear reactor core to the secondary fluid, which 
then delivers the heat to the hydrogen plant's oxygen reactor. Helium, a nitrogen/helium 
combination, or a molten salt may be used as the secondary fluid. IHX is designed to 
collect heat from the main fluid at the maximum temperature achievable (the reactor output 
temperature) for use in the hydrogen generation process. Harvego (2006) studied several 
setups for an intermediate heat exchanger between a hydrogen generation cycle and an 
NGNP. The flow rates, temperature distribution across the loops, and other IHX needs were 
also suggested based on the design, and the reactor output temperature was set at 900ºC 
in earlier research by Natesan et al. (2006). 


MATERIAL PROPERTIES OF OXYGEN REACTOR 


lkeda and Kaye (2008), as well as Trevani et al. (2011), have investigated the 
thermochemical characteristics of copper oxychloride. Because copper oxychloride is not 
commercially accessible, these two investigations discovered strategies for its manufacture. 
Ikeda and Kaye (2008) used stoichiometric quantities of CuO and CuCl, in their approach. 
To prevent the presence of moisture, the procedure was carried out in a glove box filled 
with argon. Trevani et al. (2011) used a more straightforward method. Cu,OCI, was created 
in their approach using a typical horizontal tube furnace equipped with a quartz tube and 
stainless-steel gas injection connecitors. CuCl or CuCl, was placed into an alumina crucible 
and subjected for 48 hours to a flow of high-quality dry air at 370ºC. The X-Ray Diffraction 
(XRD) analysis indicated that this approach produced very pure Cu,OCI, samples. The 
approach is more scalable, and itis capable of producing the quantities necessary for large- 
scale oxygen generation tests. Cu,OCI, is a solid that decomposes in the presence of water, 
hence it is essential to preserve it under nitrogen or dry air. 

Parry (2008) reported calorimetric measurements of copper oxychloride's heat 
capacity in the range of 20-410ºC. The Debbiee-Einstein technique was used to extrapolate 
thermodynamic data to 550ºC. Uncertainty developed in the low-temperature range due 
to the quick change in heat capacity in this area. The slope of the calorimetric specific 
heat curve drops to less than 5 J/mol.K every 100ºC above 400ºC. Differences in observed 
data were attributed to undetected heat leakage. The standard deviation of results acquired 
using this approach was calculated to be 1.1%. 

Abdulrahman (2016a, 2019a, 2020a) has derived a linear function of CuCI surface 
tension with temperature as follows: 
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0,.=0.115-5.312x 105T (4) 


Cucl 


The derivation was based on the approximate linear reduction of the surface tension 
with temperature and the equality of the surface tension to zero at the critical temperature 
of the liquid. The critical temperature that is used in deriving Eq. (4) is equal to 2435 K. 
The surface tension of molten CuCl in contact with air is equal to 0.092 N/m at the melting 
temperature of CuCl (430ºC) (Janz and Tomkins, 1983). 

There are various obstacles and challenges in utilizing the real materials from 
the thermolysis reactor products in the experimental research for the scale up analysis 
(ie. molten salt CuCl and oxygen gas). Abdulranman (2016a, 2019a, 2020a) developed 
alternative materials by simulating the hydrodynamic and heat transport characteristics 
of actual materials utilizing dimensional analyses of the Buckingham pi theorem. These 
alternate materials have been identified as liquid water at a temperature of 22ºC and helium 
gas at a temperature of 90ºC. The physical properties for He gas and liquid water are 
shown in Table 1 (Abdulranman, 2016a). Alternative materials offer a safe environment for 
experimental runs and operate at a lower temperature. Additionally, these materials are 
readily available and inexpensive. Abdulrahman's material simulation is a unique tool for 
assessing hydrodynamic and heat transport behaviours in a simulated environment prior to 
adaptation to the Cu-ClI cycle thermolysis reactor. 


Physical Property Unit Gas Phase (He) Liquid Phase (H,0) 
Temperature (T) *C 90 22 
Temperature (T) K 363 295 
Density (p) kg/m? 0.1344 998.2 
Specific Heat (C,) Jikg K 5193 4182 
Thermal Conductivity (k) Wim K 0.1687 0.6 
Dynamic Viscosity (4) kg/m s 2.267E-5 0.000975 
Molecular Weight kg/ kmol 4.0026 18.0152 
Standard State Enthalpy J/ kg mol - -2.858e+8 
Surface Tension (0) N/m - 0.0724 


Table 1: Material properties for the Helium gas and liquid Water (Abdulranman, 2016a). 


OXYGEN PRODUCTION CHEMICAL REACTION 


Experiments were carried out to determine the volatility of CuCl and the products 
of Cu,OCI, decomposition. Copper oxychloride retained its chemical characteristics 
at temperatures up to 370ºC. Cl, was released between 375ºC and 470ºC during the 
production of CuCl. Cl, then reacted with CuO to form oxygen and CuCl, as well as trace 
quantities of unreacted CuCl,. Between 475ºC and 550ºC, XRD analysis demonstrated an 
increase in the quantity of CuCI, in the products, whereas the amount of CuO decreased. 
The decomposition reaction (called path |) has the following form: 
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Cu,OCI, » CUCI + CuO + % CI, (5) 
CuO + CI, > CUCI, + 4 O, (6) 


Parallel to the competitive reaction, experimental data suggested that the secondary 
reaction was responsible for the partial drop in CuO caused by the Cl, generated in the former 
reaction. This collection was named “reaction “path |.” Another way to generate oxygen is 
by the decomposition of a combination of CuO and CuCl,, as Serban et al. investigated 
(Serban et. al., 2004). The studies were conducted in a vertical reactor at 500ºC using an 
equimolar combination of CuO and CuCl,. The XRD measurement of the solid products 
revealed a pure CuCl solid phase. The investigation used a mechanistic approach and 
determined that oxygen evolution is constrained by CuCI, decomposition. As specified in the 
following equations “path Il,” the reaction proceeded in two steps. 


CuCI, (s) > CuCl (s) + Y2 CI, (9) (7) 
CuoO (s) + Y CI2 (9) > CuCl (s) + 4% O, (9) (8) 


The competing reaction in Eg. (7) is cupric chloride's heterogeneous and endothermic 
thermolysis. Even while only 60% decomposition of CuCI, occurs at 530ºC, the addition of 
CuO seems to enhance decomposition. Additionally, it was found that the Cl, produced is 
completely consumed in the reaction with CuO, as indicated by Eq. (8), since no free Cl, was 
found in the gaseous products. The Gibbs free energy of the reaction suggests that it enters 
equilibrium around 500ºC and then proceeds spontaneously, with increasing amounts of 
chlorine available to react with CuO as the temperature rises. However, experimental data 
suggested that CuCl, does not completely decompose until it reaches 600ºC. 

Lewis et al. (2005) utilized 40 mg of laboratory-produced Cu,OCI, heated to 530ºC 
for the oxygen generation process in their experiment. Based on the quantity of copper 
oxychloride used in the experiment, it was found that the observed mass of oxygen was more 
than the theoretical limit. The presence of chlorine peaks was found to be below background 
values. Lewis et al. (2003, 2019) utilized a stoichiometric combination of CuO and CucCl, 
heated to 500ºC to achieve an oxygen output of 85%. The X-ray diffraction (XRD) analysis 
of products revealed the presence of CuCl and negligible levels of CuCI,. At500ºC, 0.02 mol 
of an equimolar combination of CuO and dehydrated CuCl, was decomposed in a separate 
experiment. The oxygen production was more than what thermodynamic equilibrium studies 
expected. The amount of chlorine released was not determined. Continuous oxygen removal 
was used in the experiments. The XRD investigation of the solid products revealed the 
presence of CuCl and CuCl,. The findings indicate that the oxygen released during copper 
oxychloride decomposition is more than the predicted yield from a combination of CuCI, and 
CuO. The reason chlorine was not found is unclear, however it might have been due to the 
presence of chlorine gas or a low yield of chlorine below the chlorine analyser's detection 
threshold. 
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Marin (2012) established a novel experimental and theoretical basis for scaling 
up a CuO*CuCl, decomposition reactor while considering the effect on the yield of the 
thermochemical copper-chlorine cycle for hydrogen production. He utilzed a Stefan 
boundary condition in conjunction with a new particle model to monitor the location of the 
moving solid-liquid interface as the solid particle decomposes due to surface heat transfer. 
Thermo gravimetric Analysis (TGA) microbalance and laboratory scale batch reactor studies 
were used to investigate the conversion of CuO*CuCI, and estimate the rate of endothermic 
reaction. At high temperatures and low Reynolds numbers, a second particle model finds 
characteristics that affect the transient chemical decomposition of solid particles embedded 
in a bulk fluid composed of molten and gaseous phases. For a particle suddenly submerged 
in a viscous continuum, the mass, energy, momentum, and chemical reaction equations 
were solved. Numerical solutions were generated and evaluated using experimental data 
on the chemical decomposition of copper oxychloride (CuO*CucCI,). 


FACTORS AFFECTING OXYGEN REACTOR SIZE 


There are many factors that affect oxygen reactor size. These factors include the 
residence time of the thermal decomposition for solid particles, terminal settling velocity of 
solid particles, production rate of hydrogen (mass balance), and reactor heating rate (heat 
balance) (Abdulranman, 2016a, 2013a). 


Terminal Settling Velocity of Particles 


The amount of time it takes for solid particles to settle in a reactor with a limited 
volume is an essential metric to consider while designing the reactor (Felice and Kehlenbeck, 
2000). The rate at which particles descend has a considerable impact on the height of the 
reactor. Itis advised that the particle fall speed be kept as low as possible in order to provide 
adequate time for the decomposition process to be completed before reaching the bottom of 
the reactor. It is not necessary to have a high reactor in this situation. 

Due to the fact that the density of the solid Cu,OCI, particle is larger than that of 
molten CuClI, the falling velocity of a single particle under gravity will be greater than that 
of molten CuCl. The particle will continue to accelerate for a brief period and then descend 
at a constant velocity, where the effective weight of the particle and the drag force are 
identical. Thus, the most essential period is the period of constant-velocity descent (terminal 
settling velocity). The terminal velocity of solid particles is determined by a number of 
characteristics, including their size, shape (roundness), solid concentration and density, as 
well as the liquid”s viscosity and density (Nguyen and Schulze, 2003). 

Particle size: In industry, solid particles are often not uniform in size, but rather 
contain a variation of sizes. It is commonly established that bigger particles settle at a 
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faster rate than smaller ones. Baldi and Alaria (1978) indicated that when dealing with a 
distribution of particle sizes, it is appropriate to describe the particle diameter as the mass- 
mean diameter when computing the settling velocity. 

Particle shape: In most theoretical computations, the form of solid particles is 
assumed to be spherical. The size of non-spherical particles may be described in terms of 
an equivalent volume of a sphere. A spherical particle's falling velocity is greater than that 
of a non-spherical particle with the same volume and density (Nguyen and Schulze, 2003). 

Particle concentration: Whena cloud of solid particles falls through a quiescentliquid, 
extra impeding actions have an impact on its falling velocity. These include the increased 
drag created by the close proximity of the particles inside the cloud and the liquid up flow 
induced by the falling particles. The impeding effects are highly dependent on the volumetric 
solid concentration in the cloud. The velocity of a hindered falling object is typically a fraction 
of the velocity of a free-falling object. For instance, with a solid concentration of 30% and 
a solid particle diameter range of 0.074-2 mm, the sand's impeded velocity is around 20%- 
40% of the single particle terminal velocity (Richardson and Zaki, 1954). 

It is important to define in the design of the oxygen reactor's size if the settling 
velocity of the solid particles allows for sufficient residence time for complete decomposition 
of these particles before they settle to the reactor's bottom. It is widely established that the 
maximum settling velocity occurs when a single spherical particle falls infinitely deep into a 
fluid. Other factors such as particle shape, reactor wall, and impeding effects all contribute 
to the particles' settling velocity being reduced. The up-flow motion of the oxygen bubbles 
also reduces the fall of particles in the oxygen reactor. 


Mass Production Rate of Hydrogen 


In the design of any chemical reactor, it is necessary to study material balances 
(mass balances) which are based on the fundamental law of mass conservation. It is difficult 
to construct or run an oxygen reactor (or any other chemical plant) in a safe and cost- 
effective manner without proper material balances on hand. The general equation of the 
material balance is (Abdulrahman, 2016a, 2013a): 


dm 


qo Mu Mo+6-—C (9) 


where: “M is the rate of change of the material m, m; and mo are the inlet and 

dt 
outlet mass flow rates respectively, G and C are the rate of generation and consumption 
respectively. Abdulrahman et al. (2016a, 2013a) investigated the flow material balance 
for continuous processes running in the steady state. To produce hydrogen, they have 


investigated the flow rates and chemical compositions of all streams entering and exiting 
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each individual piece of equipment in the Cu-Cl cycle. It has been shown that the flow rate 
of each component of the Cu-Cl cycle (in the event of full conversion of the components) 
changes linearly with the rate of hydrogen production throughout the cycle. 

Using varying residence times and hydrogen production rates, it was possible to 
compute the optimal size of an oxygen reactor. It was discovered that the size of the oxygen 
reactor changes linearly with the rate of hydrogen production and the residence time. 
According to the findings of Abdulranman et al. (2016a, 2013a), the minimal size of the 
oxygen reactor can be determined based on the mass balance of the Cu-Cl cycle, and 
the size is strongly influenced by the residence time and hydrogen production rate of the 
reactor. Also, the size of the oxygen reactor was anticipated to be large because of the 
presence of copper (Cu), which has the maximum molecular weight of 63.54 g/mol in the 
Cu-Cl cycle, and hence the largest molecular weight in the Cu-Cl cycle. 


METHOD OF HEATING OXYGEN REACTOR 

Asuitable method of heating the oxygen reactor is neededto provide enough heat that 
is necessary for solid decomposition. Heating the Cu,OCI(s) particles only is undesirable, 
because of relatively slow rate of heat transfer. A more practical and efficient option is to 
heat the molten salt inside the oxygen reactor to in turn transfer the heat within the reactor 
from the liquid (molten CuCl) to the solid Cu,OCI, (reactant) particles. The molten salt bath 
can be sustained by the reaction product itself. This approach is the most applicable and 
recommended one. 

Due to buoyancy of the gas in the molten salt, the oxygen product will leave the 
reactant particles immediately. This fast separation aids to minimize heat transfer resistance 
to the reactant particles, which then helps make the overall reaction rate closer to the 
intrinsic reaction rate. The design of the reactor requires a high efficiency of heat exchange 
and separation of reactant (solid particles) from products, as well as one product (oxygen) 
from another (molten salt). 

There are two main heating methods for reactor: electrical heating and fluid heating. 
The main disadvantage of electrical heating is the limited efficiency of electricity generation 
from nuclear heat (currently less than 35%). A more suitable method is to heat the molten 
salt by using a heating fluid (such as helium gas or molten salt) without using electricity. 

Two main configurations of molten salt heating by a fluid can be used: direct and 
indirect heating. Direct heating places the heating fluid in direct contact with the heated 
medium. The advantage of heating directly is that, heating fluid is nearly 100% efficient 
with this method. This is because all heat that is generated is absorbed directly by the 
process. This helps to speed heat-up and eliminate thermal lag. There is no intermediate 
heat transfer medium that could result in heat losses. Direct heating can be applied by using 
a Slurry Bubble Column Reactor (SBCR), where the heating gas can be introduced to the 
reactor from the bottom through a sparger. 
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Indirect reactor heating uses a heat transfer medium to deliver the heat to the reactor 
vessel. Indirect methods can be external heating of the reactor using the reactor wall (or 
immersed tubes walls). There are various advantages of indirect heating. The biggest 
advantage is that the heating fluid can typically be serviced without draining the reactor. 
Second, indirect heating often allows watt density exposed to the process fluid, which is 
the number of watts concentrated per unit surface area of process side, to be lowered by 
spreading the heat over a larger surface. Finally, overheat conditions can be limited in many 
instances by simply limiting the temperature of the heat transfer medium, so that the process 
fluid is never exposed to temperatures higher than the heat transfer fluid temperature. 

In spite ofthe above advantages, the disadvantages of indirect heating may be critical 
to the oxygen production process. The primary disadvantage is the thermal lag caused by 
using a heat transfer medium to carry the heat such as the reactor wall. The delay is caused 
by the fact that the heating fluid must first heat the heat transfer medium before the heat 
transfer medium can heat the process. If there is a large mass of heat transfer medium, 
larger heating capacities will be required to raise the temperatures. 

Different types of indirect heat transfer surfaces can be used in an industrial oxygen 
reactor such as jackets, helical tubes immersed directly inside the reactor and vertical tube 
baffles (Oldshue, 1983). A reactor jacket is usually adequate to provide the required heat 
transfer surface for low and moderate heat duties (in terms of heat duty per unit of vessel 
volume). As heat duty increases, internal heat transfer surfaces (helical coils, baffle pipes, 
or plate coils) may be required. 


Indirect Heat Transfer of the Oxygen Reactor 


Thermal scaling up of the three-phase oxygen generation reactor in the Cu-Cl cycle 
was investigated using indirect heat transfer techniques (Abdulrahman, 2016a, 2016b, 
2016c, 2019b, 2019c, 2013b, 20224). The size of the oxygen reactor necessary to generate 
sufficient heat input was investigated in these researches for various hydrogen generation 
rates. The researches employed a continuous stirred tank reactor (CSTR) that was heated 
by a half pipe jacket (Abdulranman, 2016a, 2016b) or an internal helical tube (Abdulrahman, 
2016a, 2016c). The thermal resistance of each part of the reactor system was investigated in 
order to determine its influence on the reactor's heat balance. The system thermal resistance 
was found to be dominated by the reactor wall in jacketed reactors (Abdulranman, 2016a, 
2016b), and by the helical tube wall and service side in oxygen reactor systems with an 
internal helical tube (Abdulrahman, 2016a, 2016c). 

The Cu-Cl cycle was assumed to be powered by a nuclear reactor in the above 
studies, and two different kinds of nuclear reactors were considered as a heat source for the 
oxygen reactor. The CANDU Super Critical Water Reactor (CANDU-SCWR) and the High 
Temperature Gas Reactor (HTGR) were two of these designs. It was established that a 


Engineering in perspective: Science, technology and innovation Chapter 4 


54 


larger heat transfer rate is required for CANDU-SCWR by three to four times that of HTGR, 
due to the HTGR's higher exit temperature, which results in greater temperature differential 
between the service and process sides than the SCWR. 

From a thermal balancing standpoint, the impacts of reactor diameter and aspect 
ratio on the size of the oxygen reactor were calculated. It was determined that the size of 
the oxygen reactor reduces nonlinearly when the reactor diameter or aspect ratio increases, 
with the reactor diameter having a greater influence on the reactor size than the aspect 
ratio. The pace at which the oxygen reactor's size decreases was discovered to be slowed 
by increasing the reactor's diameter or aspect ratio. 

On the service side, several kinds of working fluids were employed, including helium 
gas and CuCI molten salt, and their properties were compared. There is no discernible 
variation in the size of the jacketed oxygen reactor whether helium gas or CuCl molten 
salt is used (Abdulrahman, 2016a, 2016b). Due to the decrease in the number of oxygen 
reactors, it was established that utilizing CuCI molten salt in the service side is thermally 
superior to using helium gas in the oxygen reactor using an internal helical tube for heating 
(Abdulrahman, 2016a, 2016c). In general, both helium gas and CuCI molten salt have 
technical concerns, and both should be preserved and investigated until these technical 
difficulties are answered and an optimum decision can be chosen. 

The findings of heat transfer studies conducted on the oxygen reactor with an internal 
helical tube were compared to the results of material balances and heat transfer analyses 
conducted on the oxygen reactor with a half pipe jacket. Due to the increased size needed 
for the heat balance, it was discovered that the size of the oxygen reactor is dictated by 
the heat balance rather than the material balance. Additionally, it was shown that delivering 
heat to the oxygen reactor through an internal helical tube is more efficient than utilizing a 
half pipe jacket. 


Direct Contact Heat Transfer of the Oxygen Reactor 


In examining the direct contact heat transfer within the oxygen reactor, comprehensive 
investigations were conducted utilizing both experimental and computational fluid dynamics 
(CFD) approaches (Abdulrahman, 2015, 2016a, 2016d, 2016e, 2016f, 20169, 2018, 2020b,, 
2020c, 2022b, 2022c, 2022d, 2023a, 2023b, 2023c). These studies aimed to understand 
and analyse the heat transfer mechanisms occurring within the reactor, where oxygen gas 
interacts directly with the other components, such as molten salt or solid particles. 

The experimental studies involved setting up physical models of the oxygen reactor 
and conducting tests under controlled conditions to observe and measure the heat transfer 
rates, temperature profiles, and other relevant parameters. These experiments provided 
valuable data on how heat is transferred in the reactor, the influence of various operating 
conditions, and the effects of different reactor designs on the efficiency of heat transfer. 
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Concurrently, CFD studies were carried out to simulate the heat transfer processes 
within the oxygen reactor. Using Ansys Fluent software, detailed models of the reactor were 
created, incorporating the physical and chemical properties of the materials involved. These 
simulations allowed for a deeper analysis of the heat transfer phenomena, including the 
visualization of temperature distributions, gas holdup and flow patterns. 

By combining the insights gained from both experimental and CFD studies, 
researchers aimed to optimize the design and operation of the oxygen reactor for 
improved heat transfer performance. This holistic approach ensures a more accurate and 
comprehensive understanding of the direct contact heat transfer mechanisms in the oxygen 
reactor, which is crucial for the development of more efficient and effective thermochemical 
processes for oxygen production. 


Experimental Studies 


Scaling up the thermolysis reactor requires a thorough understanding of hydrodynamics 
and heat transfer. However, using the actual products of the thermolysis reactor, such as molten 
salt CuCl and oxygen gas, presents several challenges in experimental studies on scaling 
up. These challenges include: 1) the high melting temperature of cuprous chloride (CuCl) at 
430ºC, 2) the non-transparent dark grey colour of molten CuCl, which obstrucis the visibility 
of oxygen bubbles, 3) the highly corrosive nature of molten cuprous chloride, 4) the strong 
oxidizing property of oxygen gas, which can lead to rapid combustion of materials, and 5) the 
high-temperature process involved. To overcome these challenges, alternative materials have 
been identified using dimensional analysis to simulate the hydrodynamic and heat transfer 
behaviours of the actual materials. Liquid water at 22+2ºC and helium gas at 90+2ºC have 
been found to be suitable substitutes. These alternative materials not only provide a safer 
environment for experiments but also allow for lower operating temperatures. Moreover, they 
are readily available and more cost-effective (Abdulrahman, 2019a, 2020a). 

The experimental studies of the oxygen reactor have shed light on key aspects of 
SBCRs, including transition velocity, overall gas holdup, and direct contact heat transfer, 
providing valuable insights for optimizing their performance (Abdulrahman, 2015, 2016a, 
2016d, 2016e, 2018). In the study of the transition velocity, which marks the shift from 
homogeneous to churn turbulent flow regimes within the reactor, the experiments revealed 
that the transition velocity decreases with an increase in static liquid height and solid 
concentration. This finding suggests that the physical dimensions of the reactor and the 
density of the slurry play significant roles in determining the flow regime. Also, the research 
indicated that slug flow, a regime characterized by large gas bubbles, is not typically present 
in industrial SBCRs (Abdulranman, 2016). 

Another important aspect investigated experimentally is the overall gas holdup, 
which refers to the volume fraction of gas within the reactor. This parameter is crucial for 


Engineering in perspective: Science, technology and innovation Chapter 4 


56 


assessing the reactor's efficiency and capacity for mass transfer. The experimental study 
found that the overall gas holdup increases with the superficial gas velocity, particularly 
at lower velocities. However, it decreases with increasing static liquid height and solid 
concentration. Interestingly, the impact of solid particle diameter on gas holdup was found to 
be negligible, suggesting that particle size might not be a critical factor in certain operating 
conditions (Abdulrahman, 2016d). 

Direct contact heat transfer in SBCRs, is a process vital for maintaining optimal 
reaction temperatures and energy efficiency. The experimental results showed that both 
the volumetric heat transfer coefficient and the slurry temperature increase with the 
superficial gas velocity, with a more pronounced effect at lower velocities. Conversely, these 
parameters decrease with an increase in static liquid height and solid concentration. The 
findings imply that manipulating the gas flow rate and reactor dimensions can effectively 
control heat transfer within the system (Abdulrahman, 2015). The schematic diagram of the 
SBC experimental setup is illustrated in Fig.1 and the dimensions of the reactor column and 
sparger used in the experiments are shown in Figs. 2 and 3 respectively. 

The experimental studies provide a deeper understanding of the complex dynamics 
within slurry bubble column reactors. By elucidating the relationships between key operating 
parameters and reactor performance, these findings offer valuable guidance for designing 
and optimizing SBCRs in various industrial applications. For instance, adjusting the static 
liquid height and solid concentration can help regulate the transition velocity and gas holdup, 
while controlling the superficial gas velocity can optimize heat transfer efficiency. 
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Fig. 1 Schematic diagram of the experimental setup (Abdulrahman, 2016). 
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Fig. 2 Reactor column dimensions (Abdulrahman, 201 6a). 
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Fig. 3 Dimensions of the sparger (Abdulranhman, 2016a). 


CFD Studies 


The objective of the Computational Fluid Dynamics (CFD) simulations of direct 
contact heat transfer in the multiphase oxygen reactor, is to enhance our understanding of 
the hydrodynamics and heat transfer phenomena within these reactors. A comprehensive 
review of Eulerian CFD approaches for analysing BCR/SBCRs has highlighted the accuracy 
of these models in predicting reactor performance. Key findings include the increase in 
gas holdup with superficial gas velocity, column pressure, and gas phase density, and the 
uneven distribution of gas holdup within the reactor (Abdulranman, 2022d). 

One area of focus in CFD studies, is the volumetric heat transfer coefficient in 
SBCRs, where helium gas is injected through a slurry of water and alumina solid particles. 
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Researchers have investigated the effects of superficial gas velocity, static liquid height, 
and solid particle concentration on the volumetric heat transfer coefficient. CFD simulations 
revealed that the heat transfer coefficient increases with superficial gas velocity and 
decreases with higher static liquid height or solid concentration. These findings are crucial 
for optimizing reactor design and operation for efficient heat transfer (Abdulrahman, 201 6f). 

Temperature distribution within the slurry is another aspect that has been explored 
through 2D CFD simulations. Similar to the heat transfer coefficient, the average slurry 
temperature rises with increased superficial gas velocity and drops with higher static liquid 
height or solid concentration. However, the impact of solid concentration on temperature 
is found to be negligible. These insights are valuable for maintaining optimal reactor 
temperatures for various chemical processes (Abdulrahman, 20169, 2020b, 2020c, 2022b). 

Finally, the overall gas holdup in helium-water bubble columns has also been 
predicted using 2D and 3D CFD simulations. Results indicate that gas holdup increases with 
superficial gas velocity and decreases with higher static liquid height. The three-dimensional 
model was found to be more accurate for gas holdup in comparison to the two-dimensional 
model. This information is essential for scaling up bubble column reactors, as gas holdup 
is a key parameter describing reactor performance (Nassar, 2023; Abdulranman & Nassar, 
2023a, 2023b, 2023c). 

In Table 2, the equations utilized in the CFD analysis are displayed. Table 2 contains 
equations that are explicitly expressed in gas phase. Since the equations for the liquid phase 
are comparable to those for the gas phase, they are not repeated. Table 3 summarizes the 
setup used in the simulations. 
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Table 2: Details of equations used in the 3D CFD simulations (Abdulrahman & Nassar, 2023a). 
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Solver Type Pressure-Based 
Velocity 
General Formulation RnsoIo: 
Time Steady 
Gravity ON 
Multiphase-Eulerian 
Models Energy-On 
Viscous-Standard K - «, Standard Wall Function, Dispersed 
Water liquid 
Materials 
Helium gas 
Primary phase=liquid phase 
Phases 


Secondary Phase=gas phase 


Bubble Diameter | Sauter-mean diameter 


Scheme Phase-Coupled SIMPLE 
Gradient Least Squares Cell Based 
Momentum Second Order Upwind 

. Volume Fraction First Order Upwind 
ic Spatial Turbulent Kinetic Energy Second Order Upwind 

Discretization 
Turbulent Dissipation Rate Second Order Upwind 
Energy Second Order Upwind 


Interfacial Area Concentration | | Second Order Upwind 


Table 2 Summary of the SBCR problem setup (Abdulranman, 201 6f). 


In conclusion, CFD simulations have provided valuable insights into the hydrodynamics 
and heat transfer in slurry bubble column reactors. These advancements are instrumental in 
optimizing reactor designs and operations for various industrial applications, particularly in 
hydrogen production, contributing to the development of sustainable energy technologies. 
The results of both the experimental and CFD works are shown in Table 4. 
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Volumetric 
Heat Transfer 
Coefficient 


Average Slurry 


Temperature Overall Gas Holdup Transition Velocity 


In the industrial SBCRs, the 
slug flow regime does not 
exist and there is only one 
transition velocity, which is 
between the homogeneous 
and heterogeneous flow 
regimes. 


Increases by increasing the superficial gas velocity with a higher rate 
of increase at lower superficial gas velocity. 


Decreases by increasing the static liquid height. 


Decreases by increasing 
solid concentration 

and at a higher solid 
concentration, the 
changing rate of the 


The decrease of the overall gas holdup 


paia d OR with the superficial gas Decreases by increasing 
concentration. concentration is velocity and/or the solid solid concentration. 
negligible. concentration is lower. 
The effect of the solid 
particle diameter on 
overall gas holdup is 
negligible. 
The distribution along 
the cross-section of 
The rate of decrease with the solid the column is unequal, 
concentration are approximately the same | where the gas holdup is 
for different superficial gas velocities. higher at the centre of the 


column and lower near 
the wall region. 


Profiles calculated from CFD models, generally under-predicted the 
experimental data of Abdulrahman (2015, 2016d, 20161). 


The CFD models correctly predicted the effects of superficial gas 
velocity, static liquid height and solid concentration. 


The CFD results were validated for superficial gas velocities up to 0.15 
m/s, and aspect ratios up to 4. 


Table 4: Results of experimental and CFD analyses for direct contact heat transfer (Abdulrahman, 
20220). 


CONCLUSIONS 


Recent improvements in the multiphase oxygen reactor for the copper-chlorine 
thermochemical cycle were presented and assessed in this review work. The study 
identified previous experimental and theoretical findings acquired from research groups 
working on the development of scaled-up oxygen reactor facilities. The oxygen reactor's 
type and description were discussed in details and it is concluded that using SBCR is the 
best option for the oxygen reactor. Oxygen production chemical reaction was described in 
details and it was found that the reaction is an endothermic reaction with a reaction heat 
of 129.2 kJ/mol and a temperature of 530ºC which is the highest temperature in the Cu-Cl 
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cycle. The source of heat can be from nuclear reactors or solar energy sources. Also, the 
material properties of the oxygen reactor contents were extensively explained. Moreover, 
the factors that affect the size of the oxygen reactor were discussed. The factors include the 
residence time of the thermal decomposition for solid particles, terminal settling velocity of 
solid particles, production rate of hydrogen, and reactor heating rate. The techniques of heat 
transfer required by the oxygen reactor, such as indirect and direct contact heat transfer, 
were also analysed in this study. It was concluded that using direct contact heat transfer is 
preferred to indirect heat transfer methods. 
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